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In  this  work,  we  consider  a  stock  price  process  subjected  to  idiosyncratic  Levy  jumps  and 
global  structural  changes  attributed  to  interventions  due  to  a  semi-Markov  process.  The 
semi-Markov  process  decomposes  both  the  time  and  state  domains  of  the  price  process 
into  sub-intervals  and  price  state  sub-domains  respectively,  where  a  Levy-Ito  process 
operates.  The  Levy  jumps  decompose  the  space  domain  of  the  currently  operating  Levy 
process.  We  derive  an  infinitesimal  generator  for  a  stock  price  process  and  a  closed  form 
expression  for  the  conditional  characteristic  function  of  a  log  price.  The  former  result 
is  used  to  derive  a  PIDE  satisfied  by  option  prices,  while  the  latter  could  be  used  to 
retrieve  risk  neutral  densities  via  Fourier  transform  and  price  European  vanilla  options. 
In  the  sequel,  we  derive  the  characteristic  function  of  the  residence  time  of  a  semi-Markov 
process.  Incompleteness  of  the  market  is  exhibited  through  a  general  change  of  measure. 
For  pricing  purpose,  the  minimum  entropy  martingale  measure  is  defined  as  an  Esscher 
transform. 

Keywords:  Semi-Markov  process;  regime  switching  models;  characteristic  function; 
calibration  and  simulation  option  prices;  minimum  entropy. 


1.  Introduction 

The  well-known  model  developed  by  Black  &  Scholes  (1973),  despite  its  slew  of 
laureates  has  long  shown  well  documented  weaknesses  in  keeping  up  with  the  styl¬ 
ized  facts  of  financial  asset  returns  and  option  prices.  Smiles,  smirks  and  skews 
are  well  documented  empirical  features  of  implied  volatilities  (Bulla  2006,  Tankov 
2003,  Elliott  et  al.  2005,  Jackson  &  Jaimungal  2009),  unexplained  in  the  context 
of  Black-Scholes  model.  Moreover,  stylized  facts  of  financial  time  series  also  cast 
a  doubt  on  the  appropriateness  of  the  normal  log  return  distribution  assumption. 
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The  literature  goes  about  solving  these  issues  in  two  main  ways.  The  first  approach 
uses  time  dependent  deterministic  volatility  models  (Kruse  2003,  Han  2006)  to  cap¬ 
ture  most  of  option  market  empirical  properties.  However,  it  has  been  shown  by 
Gonzalez-Urteaga  (2012)  that  risk  neutral  volatilities  behave  in  a  random  man¬ 
ner.  This  leads  to  the  development  of  the  second  modeling  approach  consisting  of 
stochastic  volatility,  local  volatility  and  regime  switching  models.  Stochastic  volatil¬ 
ity  models  (Hagan  et  al.  2002,  Heston  1993)  are  based  on  the  assumption  that 
volatility  is  a  dynamic  process.  In  local  volatility  models  (Beckers  1980,  Dupire 
1997),  the  volatility  depends  on  time  and  stock  price  through  a  deterministic  func¬ 
tional.  In  both  cases,  in  addition  to  the  possibility  of  a  misspecified  functional  form 
of  the  volatility  (Chourdakis  2005),  the  volatility  surface  often  lacks  smoothness  and 
at  times,  takes  nonsensical  and  counterintuitive  forms.  One  of  the  main  advantages 
of  regime  switching  models  as  noted  by  Chourdakis  (2005),  is  the  interpretability 
of  the  market  states  while  disassociating  with  the  very  restrictive  functional  form 
assumption  of  the  local  and  stochastic  volatility  models. 

The  present  paper  is  an  attempt  to  extend  the  latest  semi-Markov  switching 
models  for  stock  price  and  at  setting  up  a  general  theoretical  framework  to  study 
qualitative  and  quantitative  properties  of  asset  price  processes.  To  the  best  of  our 
knowledge,  mostly  Markovian  interventions  on  price  processes  are  studied  in  Chan  & 
Zhu  (2014b, a),  Hainaut  &  Le  Courtois  (2014),  Hainaut  (2011),  Momeya  (2012), 
Jackson  &  Jaimungal  (2009)  and  Chourdakis  (2005,  2002).  Recently,  stochastic 
models  (Preda  et  al.  2014,  Swishchuk  &  Islam  2011,  Hunt  &  Devolder  2011,  Hunt  & 
Hahn  2010,  Ghosh  &  Goswami  2009)  under  the  influence  of  a  semi-Markov  process 
have  been  examined.  A  Stochastic  Maximum  Principle  for  semi-Markov  switching 
jump  diffusion  models  (Deshpande  2014)  has  been  established,  leaving  out  the  class 
of  Levy  processes  with  infinite  activity. 

The  organizational  outline  of  the  paper  is  as  follows:  In  Sec.  2,  we  introduce 
the  necessary  definitions  and  notations  and  we  present  known  results.  In  Sec.  3,  we 
find  a  closed  form  solution  of  a  Levy  type  of  SDE.  In  Sec.  4,  we  derive  Ito  differ¬ 
ential  formula  and  the  infinitesimal  generator  for  a  class  of  stochastic  linear  hybrid 
models  under  semi-Mar kovian  and  Levy- type  structural  perturbations.  Section  5 
is  concerned  with  the  derivation  of  a  closed  form  characteristic  function  of  the  log 
price  process.  This  is  useful  for  recovering  risk  neutral  densities  to  estimate  option 
prices.  Moreover,  this  provides  an  alternative  tool  to  the  computationally  extensive 
continuous  time  MCMC  and  the  two-step  numerical  integration  procedure  (Hunt  & 
Hahn  2010,  Ghosh  &  Goswami  2009),  for  simulating  option  prices  and  calibrating 
model  parameters  to  market  option  prices.  In  Sec.  6,  we  exhibit  a  general  change  of 
measure  and  two  risk  neutral  measures  of  interest,  namely,  the  minimum  entropy 
martingale  measures  and  the  pricing  kernel  discussed  in  Siu  &  Yang  (2009).  The 
latter  accounts  for  the  regime  risk,  the  jump  risk  and  the  Levy  risk.  Section  7  is 
devoted  to  the  presentation  of  a  couple  of  option  price  formulas.  The  first  is  the  well 
known  Fourier  transform  method  based  on  the  method  of  Carr  &  Madan  (1999). 
The  second  formula  is  a  slight  modification  of  the  integral  formula  developed  in 
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Ghosh  &  Goswami  (2009).  In  Sec.  8,  numerical  illustrations  are  given  to  exhibit  the 
usefulness  of  the  presented  results. 

2.  Preliminary  Definitions  and  Results 

In  this  work,  T*  <  oo  and  T  G  [0,T*].  T*  and  T  stand  for  the  time  horizon  of 
the  market  and  the  maturity  time  of  a  contingent  claim,  respectively;  (fi,F,  P)  is  a 
complete  probability  space;  0  is  a  semi-Markov  process  defined  on  M+  x  P) 

into  P,  where  P  is  an  at  most  countable  subset  of  the  set  of  natural  numbers  N  and 
M+  =  [0,  oo).  For  each  n  G  N,  Tn  stands  for  the  nth  jump  time  of  9.  For  s  G  [0,  T] 
and  0S-  =  j,  (Zd_)se[0jT]  is  the  Ito  Levy  process  with  small  and  big  jumps  G  and 
H  and  Levy  triplet  (/i(j),  cr(j),  i/(j, )),  where  /x(j),  a(j)  and  n(j, )  are  the  drift  rate, 
the  diffusion  rate  and  the  Levy  measure,  respectively.  (ht)te[o,T]i  (^t)te[o,T]  and 
(Hit  V  Lt)te[o ,t]  are  sub-algebras  of  F  generated  by  the  collection  of  Levy  processes 
PJ,  V  j  G  P,  the  semi-Markov  process  0  and  (T^ ,  0t)te[o,T],  respectively.  Let  (/?n)n> o 
and  (®t)te[o,T]  be  a  discrete-time  real  valued  stochastic  process  and  the  sub-sigma 
algebra  of  F  adapted  to  the  discrete  process  /?n,  respectively.  We  denote  the  enlarged 
filtrations  (Lt  V Bt)tG[0jT]  and  (ItV® t)te[o,T]  by  (L*t)tG[o,T]  and  (fflt)te[o,T],  respec¬ 
tively.  Let  -,  •)  be  the  Poisson  random  measure  with  compensator  v(j,dz).  It 
is  also  assumed  that  the  sequence  (/3n)n> o  is  independent  of  both  •,  •)  and  the 
Brownian  process  Bt ,  for  j  G  P. 

Definition  2.1.  (Cinlar  1969)  Let  6  and  {Tn}^=1  be  a  semi-Markov  process  and 
its  jump  time  sequence  with  To  =  0,  respectively.  A  couple  ( 0n,Tn )  is  called  a 
Markov  renewal  process  with  kernel  Q  induced  by  the  semi-Markov  Process  (0*),  if 
it  satisfies: 


P(9n  =  j,Tn  <t  \  (6k,Tk),k  =  1,2, . . .  ,n  -  1) 

=  P(9n  =  j,Tn  -  Tn- 1  <t  —  Tn-i  |  0n-i,Tn-i) 

=  Q(en.1,j,t-Tn-1),  (2.1) 

where  0n  stands  for  0^. 


Remark  2.1.  rn  =  Tn+i  —  Tn  denotes  a  holding  time  at  Tn.  The  holding  times 
conditional  on  the  current  state  are  independent  Cinlar  (1969).  The  kernel  in  (2.1) 
can  be  represented  as 


Q(i,  T  t-Tn)  =  P{0n  =  j,  Tn  -  Tn_!  <  t  -  Tn_!  |  Qn—\  =  i).  (2.2) 

Moreover,  for  ( 9n,Tn )  =  (0„(t),  Tn(t)),  Vf  G  [0,T],  where 

n(t)  =  max{fc  <  t}.  (2-3) 

In  particular, 

=  P(0n  =  j,Tn  -  Tra_i  <  1 1  9n- 1  =  i ).  (2.4) 
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Furthermore,  we  define 

Pij  =  Jim  Q(i,j,t),  (2.5) 

t — XX) 

where  is  called  the  steady  state  transition  probability  of  the  embedded  Markov 
chain  from  state  i  to  state  j  with  i,  j  E  E  and  n(E)  =  m. 

For  the  sake  of  completeness,  we  present  survival  distribution  and  sojourn  time 
distributions  associated  with  the  semi-Markov  process  0. 


Definition  2.2.  The  conditional  cumulative  distribution  of  the  holding/sojourn/ 
residence  time  (respectively,  survival  function)  given  that  0  transits  from  a  state 
i  to  state  j  is  defined  by  F(t\i,j)  =  P(rn  <  t\0n  =  j,6n-\  =  i)  (respectively, 


In  the  following  lemma,  we  outline  a  few  well  known  properties  of  semi-Markov 
processes  (Cinlar  1969,  Ghosh  &  Goswami  2009). 


Lemma  2.1.  The  kernel  of  the  semi-Markov  process  defined  in  (2.4)  is  represented 
by 


Moreover , 


)  =Pi,jF(t\i,j). 


S{t\i)  =  l-  ^2pijF(t\i,j), 

jeE 

Hr  |  flp  =  t)  =  -fy  Ki) 
P(Ti  >  s  I  Oo  =  i)  S(s-  I  i)  ’ 


—  Pij 


zMj  ij) 

S(t-\i ) 


for  r  >  s, 


(2.6) 

(2.7) 

(2.8) 

(2.9) 


Proof.  We  first  establish  (2.6).  From  (2.4)  we  have 

=  P(0n  =  j,Tn  -  T„_ i  <  1 1  9n- 1  =  i), 

=  P(0n  =  j  I  0n- 1  =  i)P(Tn  -  Tn- 1  <  1 1  0n- 1  =i,0n  =  j ) 


=  PijF(t\i,j). 

For  the  proof  of  (2.7),  we  use  Definition  2.2  and  (2.6). 

F(t  |  i)  =  P(rn  <t\0n=i ) 


=  PJPi,iFd\hj)- 

j£E 

Hence, 

S(t\i)  =  1  -^PijFit  \i,j). 

j£E 


(2.10) 
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This  completes  the  proof  of  (2.7).  For  the  proof  of  (2.8),  for  r  >  s,  we  have 

f(r  |  0o  =  i) 


P(Ti  >  s  |  0o  =  i)  ’ 


(definition  of  conditional  density) 


-f  (r|0 


p(Ti  >  s  |  es  =  i)  ’ 


(definition  of  survival  function) 


-f(H  o 


(2.11) 


This  completes  the  proof  of  (2.8).  The  proof  of  (2.9)  follows  from  the  definition  of 
Hazard  functions  and  (2.6).  This  completes  the  proof  of  the  lemma.  □ 


Remark  2.2.  A  homogeneous  Markov  process  is  a  particular  case  of  semi-Markov 
process.  Hence,  Qij(t)  =  Pij(l  —  eq M).  We  also  have  the  following  relationship: 

qij=Pijq(i),  (2.12) 

where  (qij)mxm  is  the  Infinitesimal  generator  (intensity  matrix)  of  a  Markov  pro¬ 
cess,  and  (pij)mxm  is  its  transition  probability  matrix  defined  in  (2.5).  The  following 
can  also  be  inferred  from  (2.6)  and  (2.7): 

5(t|*)  =  l-^Q(*,i,t).  (2.13) 

j£E 

Definition  2.3.  Let  yt  be  the  backward  recurrence  time  of  the  semi-Markov  process 
0  at  time  t.  yt  is  defined  as  follows: 

yt  =  ^(t  -  rn)i(Tn<t<Tn+1),  (2.14) 

n>  0 

where  the  sequence  {Tn}^L0  is  introduced  in  Definition  2.1. 

Definition  2.4.  Let  ^  x  M  x  R  be  the  random  Poisson  measure  with 

intensity  measure  z/,  H  and  G  smooth  functions  defined  on  M+  x  M  into  M,  with 
G  satisfying  the  condition:  fz€R((l  +  -H'2(z,.?'))l|*|>1  +  G2(z,j)l\z\<1)v(j,  dz)  <  oo, 
V  j  G  E.  Moreover,  i/j  =  ^  —  v  denotes  the  compensated  Poisson  measure  associated 
with  ip. 

In  the  following  we  present  a  lemma,  which  would  be  used,  subsequently. 

Lemma  2.2.  Let  (anibn)  and  (cn,dn)  be  two  renewal  processes  defined  on  the 
same  probability  space  (fi,F,  P )  and  state  space  E.  Then  the  renewal  processes  have 
identical  transition  probability  matrices  and  sojourn  time  distributions ,  respectively. 


Proof.  From  (2.5),  it  is  clear  that  the  transition  probability  and  the  holding  time 
distribution  are  completely  defined  by  the  kernel  matrix.  In  fact,  we  have, 

lim  Q(i,j,  t )  =  lim  PijF(t  \  i,j)  =  pijt 

t — >oo  t — >oo 
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and  hence, 

Plj 

This  establishes  the  result.  □ 

Lemma  2.3.  Let  n{t)  be  defined  as  in  (2.3).  The  pair  (Ot,t  —  Tn is  a  Markov 
process. 

Proof.  Let  be  s  <  t  with  t  G  [Tn,Tn+i)  and  s  G  [Tm,Tm+i)  (m  <  n).  For  u  <  s, 
we  have: 

P(9t  =  i,t  —  Tra(t)  <  a  |  (9U,  u  -  Tn (for  some  a  G  M+) 

=  p(^Tn  =  i,Tn+1  —Tn<a  \  ( 9u,u  -  Tn(u))),  (for  t  G  [Tn,Tn+ 1)) 

=  P(^Tn  =  i,Tn. |_i  —  Tn  <  a  |  (9rk,Tk),k  <  m),  (for  s  G  [Tm,Tm+ 1)) 

=  P {0rpri  =  i,  Tn+ 1  —  Tn  <  a  |  #Tm),  (Markov  renewal  process  property) 

=  P(#t  =  i,  t  —  Tn(£)  <  a  |  #s),  (definition  of  n(t)). 

Hence,  the  probability  at  a  future  time  depends  only  on  the  most  current  informa¬ 
tion  at  time  s.  This  shows  that  (0t,t  —  Tn(t))  is  a  Markov  process.  □ 

Remark  2.3.  For  the  remainder  of  this  paper  0  is  a  semi-Markov  process  with  jump 
time  Tn,  with  sojourn  time  rn  =  Tn+1  -  Tn  ~  /(|0n,0n+i)  with  CDF  F(|0n,  0n+i) 
and  with  survival  CDF  S(|0n ,  0n+ 1)  =  1  —  F(|0n,  0n+i)-  The  semi-Markov  kernel  is 
denoted  Q(i,  j,  t),  the  backward  recurrence  time  yt  is  defined  in  (2.14)  and  (pij)mxm 
is  the  transition  probability  matrix  of  the  embedded  Markov  chain. 

3.  Method  for  Finding  Closed  Form  Solutions 

In  this  section,  we  find  a  closed  form  solution  of  a  Levy-type  Linear  Stochastic  Dif¬ 
ferential  Equation  under  semi- Markovian  structural  perturbations.  The  presented 
extension  is  based  on  the  procedure  described  in  Ladde  &  Ladde  (2013).  The  useful¬ 
ness  of  the  result  is  at  least  two-fold.  It  is  used  to  establish  the  martingale  property 
for  certain  processes  in  Sec.  6.  In  addition,  it  is  used  to  formulate  a  general  expres¬ 
sion  for  the  simple  return  process  with  any  Levy  and  semi-Markov  jump  choices. 
We  consider  the  following  Levy- type  SDE: 

dxt  =  xt-dL°t ,  x(0)  =  xo,  (3.1) 

where 

dLet  =  n(6t)dt  +  (j{6t)dBt  +  f  H{z,6t)fi{0udz,dt) 

J\z\>l 

+  [  G(z,0t)^(Ot,dz,dt),  (3.2) 

J\z\<l 
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0  is  the  semi-Markov  process  defined  in  Sec.  2;  v  G  and  H  are  in  Definition  2.4. 
Following  the  procedure  described  in  Chapter  2,  Ladde  &  Ladde  (2013),  we  break 
down  (3.1)  into  the  following  four  types  of  simplified  SDEs: 

dx \  =  x\fi[di)dt 
dx \  =  x^-cr(6t)dBt 

<  dx\  —  x\-  /  H(z,0t)^>(0udz,dt)  (3.3) 

J\z\>l 

dx\  —  x\-  f  G(z,0t)^>(0udz,dt). 

J\z\<l 

Imitating  the  procedure  in  Ladde  &  Ladde  (2013),  the  closed  form  solution 
processes  of 

dx\  =  x\- ii{6t)dt  and  dx\  =  a(0t)x^dBt 


are 


x,  = 


xt  =  exp 


ci  and 


/  n(0s)ds 

i  a\es)ds+  f  a(6s)dBs 

*  Jo  Jo 


(3.4) 


C2, 


respectively;  ci  and  C2  are  arbitrary  constants.  We  next  consider  the  third  type  of 
SDE  in  (3.3) 


dx\=x\-  [  H(z,9t)'ip(0t,dz, 

J\z\>l 


dt). 


We  seek  a  solution  of  a  form 


xt  =  exp 


f  f  f4{z,6s)'tp(9s,dz,ds) 
Jo  J \z\>l 


C3, 


(3.5) 


(3.6) 


where  is  an  unknown  smooth  function  to  be  determined,  and  C3  is  a  real  random 
variable.  The  Ito  integral  for  pure  jump  processes  (3.6)  yields 


Lt+At 


~xt  =  E  (Xs~Xs~) 


t<s<t+At 
pt~\~At 


=  /  (* 

Jt  J\z\>l 

rt+At  r 

=  A-  /  ( 

Jt  J\z\>l 


3  Ja(z,8s-)  _  3 


■_  eu*K~i-s- /  _  x*_  )%j)(ds-  ,  dz,  ds) 


[eMz,e.-)  _  l )tp(es-,dz,ds). 

As  At  becomes  very  small,  (3.7)  reduces  to 

Mz,et)  _  1  )${pudz,dt). 


dx}  =  x /  (e 
J\Z\>1 


(3.7) 


(3.8) 
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Since  x 3  is  solution  of  stochastic  differential  equation  (3.5),  we  repeat  the  procedure 
described  in  Ladde  &  Ladde  (2013)  and  obtain 


exp(/4(2,0t))  -  1  =  H(z,0t). 

Hence, 

h(z,0t)  =  ln(l  +  H(z,0t)). 

Therefore,  the  general  solution  of  (3.5)  is  represented  by 


x3  =  exp 


[  [  ln(l  +  H(z,6s))'ip(0Sldz,ds) 
JO  J\z\>  1 


C3- 


(3.9) 

(3.10) 

(3.11) 


x3  is  almost  surely  finite.  Finally,  we  find  a  solution  of  the  following  stochastic 
differential  equation 


dx\  =  x\-  f  G(z,0t)'ip(0t,dz,dt). 
J\z\<l 

We  seek  a  solution  process  of  (3.12)  in  the  following  form: 


xt  =  exp 


[  [  fo(z,6s)*l)(6s-,dx,ds)  +  [  [  fo(z,  6S- )v(6s- ,  dz)ds 

JO  J\z\<  1  Jo  J\z\<l 


(3.12) 

C4, 

(3.13) 


where  and  /@  are  unknown  smooth  functions  to  be  determined,  and  C4  is  a  real 
valued  random  variable.  xA  in  (3.12)  is  an  exponential  function  semi-martingale  of 
the  form  V  =  IoI\z\<l  f5(z,0s)rp(ds,dx,ds)  +  /0*  fe(z,  9s)i/(0s,  dz)ds.  Applying 
the  Ito  formula  for  discontinuous  semi- martingales,  Applebaum  (2009),  we  have 

dx4  -  tp±dvc  ,  ppd  —  —  (  *  *  dxA 

dxt  -  gvdv  +  2  Qv2 


d(vc)d(vc)  +(Ax$  ~  — AL 


+  xt~ 


[  (fe(z,9t-)  +  exp(f5(z,0t-))  -  1  -  h{z,et-))v(0t-  ,dz)dt 

[  (exp(f5(z,0t-)) -l)ij>(0t-,dz,dt).  (3.14) 

J\z\<  1 


Again,  following  the  procedure  for  finding  solution  processes  in  Ladde  and  Ladde 
(2013),  we  get 


Hence, 


\f6(z,0t)  +exp (f5(z,9t))  -  1  -  fs(z,6t)  =  0 
[exp (f5(z,0t))  -  1  =  G{z,0t). 

f  fb{z,  0t)  =  ln(l  +  G(z,  0t)) 

l/eMt)  =  ln(l  +G(z,0t))  ~  G(z,6t)- 


(3.15) 


(3.16) 
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Therefore, 


xt  =  exp 


[  [  [ln(l  +G(z,Qs-)) -G(z,6s-)]v(6s-,dz)ds 
JO  J\z\<l 

+  f  f  ln(l  +  G(z,  6S-))^{6S- ,  dz,  ds) 

JO  J\z\<l 


(3.17) 


The  product  of  x 1  and  x2  from  (3.4)  analog  with  x 3  and  x 4  in  (3.11)  and  (3.17), 
respectively,  yields  the  solution  of  initial  value  problem  (3.1) 


xt  =  xq  exp 


/* 


-  ^2{0S~) 


+  [  [ln(l  +G(z,6s-))  -  G(z,6s-)]v(6s-,dz)  ds 

[  a(0s-)dBs+  [  f  ln(l  +  G(z,6s-))ip(6s-,dz,ds) 
Jo  Jo  j\z\<  1 


+ 


+ 


tr  ln(l  +  H(z,0a-)M0a-,dz,ds) 
Jo  J\z\>l 


(3.18) 


In  the  following,  we  present  a  few  versions  of  (3.18). 

Remark  3.1.  We  note  that  adding  and  subtracting 

[  [  G(z,  0S- )4>{0S- 1  dz,  ds)  and  f  f  H(z,98-)ip(68-,dz,ds), 
Jo  J\z\<l  Jo  J\z\>l 

(3.18)  reduces  to 

xt=x0ex.p  /  ii(6s-)ds+  /  cr(6s-)dBs+  /  /  G(z,0s-)^(0s-  ,dz,ds) 

Jo  Jo  Jo  J\z\<l 

+  ft  H{z,es-)^{es-,dz,ds)-  \  f  a\es-)ds 

Jo  J |d>i  ^  Jo 

+  [[  [ln(l  +  G(z,6s-))-G(z,0s-)]i'(9s-,dz)ds 
Jo  j\z\<l 

r  [  Hi  +  G(z, Qs-))~  G(z , ea- Mes- , dz) 

Jo  J\z \<1 


+ 


+ 


r  f  Hi  +  0s-))~  H(z ,  0S- )}^{0S- ,  dz , 

Jo  J\z\>l 


ds) 


(3.19) 
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Moreover,  adding  and  subtracting 

pt 


[  [  G(z,0s-)i>(63-,dz,ds), 

JO  J\z\<l 

nH(z10s-)2p(0s-  ,dz,ds),  and 
,  s|>i 

[  [  H(z,  6S- )v(Qs- ,  dz)ds, 

Jo  J\z\>\ 


(3.18)  becomes 


xt  =  xq  exp 


+ 


+ 


[  p(0s-  )ds  +  f  a(6s-)dBs+  f  [  G(z,0s-)^{0s-  ,dz,ds) 

Jo  Jo  Jo  j\z\<l 

+  [  [  H(z,0s-)$(6s-,dz,ds) 

Jo  J  |  |  >  1 

tf  HH(z,es-)  +  l)v{es-,dz)ds-1-  [*  a2(0s- )ds 
Jo  J\z\>l  Z  Jo 

r  [  Hi  +  G(z, Qs-))~  G(z , 0S- Mes- , dz, ds ) 

Jo  j\z\<l 

+  [  f  [ln(l  +  G(z,6s-))  —  G(z,6s-)\v(6s- ,dz)ds 

JO  J\z\<l 

+  [  [  I111!1  +  H(z,08-))^H (z,  0s-  )]$(68- ,  dz,  ds) 

Jo  J\z |>1 


(3.20) 


In  the  following  remark,  we  take  a  look  at  a  few  particular  cases  of  interest  which 
will  be  used,  subsequently. 

Remark  3.2.  If  H(z,6s),  G(z,6s)  and  Les  in  (3.2)  are  replaced  by  eH(z'Gs^  —  1, 

eG(z,0s)  _  ^  anq 

dLes  =  n(9.)da  +  <j{Qs)dBs  +  f  [, sG (*>«■)  -  l]$(0a,  dz,  ds) 

J\z\<  1 


7  [■ 


M<i 
_  1  ]ip(6s,dz,ds), 


(3.21) 


respectively,  then  the  solution  of  the  IVP  (3.1)  in  (3.18),  (3.19)  and  (3.20)  reduce 
to 

pt 


xt  =  xq  exp 


-  \c2(6s-) 


+ 


f 

[  [G(z,es-)  +  1  -  eG('z’e‘-)]v(6s-,dz) 

J\z\<l 


ds 
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+ 


+ 


xt  =  xq  exp 


f  a(9s-)dBs  +  f  [  G(z,9s-)ip(9s-,dz,ds) 

Jo  Jo  J\z\<l 

[  [  H(z,93-)i/j(93-,dz,ds)  , 

Jo  J\z\>l 

[  ii(93-)ds+  [  a(9s-)dBs 
Jo  Jo 


+ 


// 

Jo  j\z\<l 

ft 

Jo  J\z 


[< eG(z’ff° -)  -  1  }tp(9s-,dz,ds) 


0H(z,es-)  _ 


1  ]ip(Os-,dz,ds) J  a2(9s-)ds 


+  f  f  [G(z,8s~)  —  eG^z,9s~^  +  l]v(03- ,dz)ds 
Jo  J\z\<  1 

f  f  [G(z,93~)  -  eG(-z'e‘-'>  +  l]$(6s-,dz,ds) 
Jo  J\z\<l 

f  !  [H(z,  9„- )  -  eH^’es- >  +  1  }iP(9s- ,  dz,  ds ) 

Jo  J\z\>l 

[  fi(93-)ds  +  [  a(9s-)dBs 
Jo  Jo 


+ 


+ 


xt  =  xq  exp 


leG(z,0s-)~  ^(9s-,dZ,ds) 


[eH(z’es- 1  —  1  ]ip(Os-  ,dz,ds) 


+fl 

Jo  j\z\<  1 

+// 

Jo  J\z . 

+  f  f  H(z,9s-)v(6s-,dz)ds  -  ^  [  a2(0s-)ds 

Jo  J |z|>l  ^  Jo 

+  [  f  [G(z,93~)  -  eG(-z'e‘-'>  +  l]-ip(9s-,dz,ds) 

Jo  J\z\<l 

+  f  [  [G(z,6s~)  —  eG^z,6s~^  +  l\i/(0s- ,dz)ds 

JO  J\z\<l 

+  [  [  lH (*>  °s- )  -  )  +  1  ]$(08- ,  dz,  ds) 

Jo  J\z\>l 


(3.22) 


(3.23) 


(3.24) 
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In  addition,  if  p(0s)  in  (3.2)  is  replaced  by  [/jl(Os)  +  \<j2(0s)  +  f\z\d[eG(z,es>)  —  1  — 
G(z,  0s)]v(6s,  dz)],  then  (3.22),  (3.23)  and  (3.24),  respectively,  reduce  to 

[  p(Os-)ds+  [  cr(Os-  )dBs 
Jo  Jo 


xt  =  xq  exp 


+  [  [  G(6S- ,  z)'i/j(6s- ,  dz,  ds) 

Jo  J\z\<l 

+  [  [  H(0S- ,  z)^{6s- ,  dz,  ds) 

Jo  J\z\>l 


j\z\>  X 
=  %o  exp [Let], 


xt  =  xq  exp 


[  n(6s-)ds+  f  a(0s-)dBs 
Jo  Jo 


+ 


+ 


/'/ 

Jo  j\z 

fj 

Jo  j\z 


pGW,-)  - 


N<i 


5  ^2,  ds) 


/o  ^p|>l 

nt 


\eH(z,es-)  _  i ]i/>(ea-,dz,ds) 


xt  =  xq  exp 


f  [  [G(z,es-)-eG^-Ames-,dz) 

Jo  J |z|<l 

f  f  [. H(z ,  es- )  -  eH^s- )  +  I]^es- ,  dz,  ds) 

Jo  J |z| >1 

[  n{6s-)ds+  f  o(8s-  )dBs 
Jo  Jo 


+ 


+ 


fj 

Jo  J\z 

fj 

Jo  J \z 

fj 

Jo  J\z\ 


.GWs-)  - 


1  ]VK$s-  5  dz,  ds) 


[eH(z,es-)  _  i }^p(es-,dz,ds) 


/ 0  J\z\>l 

rt 


\eHWs-)  -i]u(es-,dz)ds 


+  [  f  [G(z,0a-)  ~eGG’es-)  +  l]t/}(es-,dz,ds) 

Jo  J\z\<l 

+  [  [  [■ H(z ,  es- )  -  eH^s- )  +  i]^es- ,  dz,  ds) 

Jo  J |z| >1 


(3.25) 


(3.26) 


where  Le  is  defined  in  (3.2).  Moreover,  if  fi(Ot- )  in  (3.1)  is  replaced  by  [l40s- )  + 
icr2(6»s-)  +  /|z|<1[eG(z,^“)  -  1  -  G{z,0s-)]u(0s-,dz)  +  Jjz|>1  H(z,  0a-)v{08-  ,dz)\, 
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then  (3.24)  reduces  to 

/i(6s-)ds+  [  a(Os-)dBs 
Jo 


+ 1 

tl 

f  [eG(*,e.-)-i]^(0a_}dZ)ds) 

J 

0  J 

rt 

M<i 

r 

+ 1 

f  I 

1  [eH(z’e«~)  -l]iP(0s-,dz,ds) 

0  J 

W>i 

+ 1 

t'j 

f  [G(z,9s~)  -  eGtzfis-)  +  l]$(6s-,dz,ds) 

J 

M<i 

+  1 

fj 

[  [H(z,  9S- )  -  eH^e>~ >  +  1$(6>,_ ,  dz,  ds ) 

J 

M>  i 

Xt  =  xq  exp 


/ 


(3.27) 


4.  Ito  Differential  Formula 


In  this  section,  we  define  the  asset  price  model  and  we  derive  the  infinitesimal 
generator  of  the  quadruplet  (t,y,0,x).  We  denote  L 9  the  Ito  Levy  process  with 
small  and  big  jumps  G  and  H  and  Levy  triplet  (/i(0*  ),<j(0t),v(0t,dz))  defined  (3.2). 
Following  the  argument  used  in  Ladde  &  Ladde  (2013),  we  define  a  linear  stochastic 
hybrid  dynamic  model  for  stock  price  process  under  structural  perturbations  of 
semi-Markov  and  Levy  processes. 


Definition  4.1.  A  linear  stochastic  hybrid  dynamic  model  under  semi-Markov  and 
Levy  structural  perturbations  is  defined  as  follows: 

(dx(t)  =  x(t~)dL9tn,  x(Tn)  =xn,  t  e  [Tn,Tn+i)  ^ 

=  (JnX(Tn  ,  Xn_i,  Xn_i),  x(0)  =  Xo,  Tl  G  7(1,  Oo)  =  N, 

where  {Tn}™=1  is  an  increasing  sequence  of  jump/regime  switching  times  of  the 
semi-Markov  process  6  with  To  =  0  introduced  in  Definition  2.1;  for  n  G  /( 0,  oo)  = 
{0, 1,2,3,...},  f3n  denotes  the  discrete  time  state  jump  process  caused  by  the  semi- 
Markov  process  from  state  0n-i  at  Tn_i  to  0n  at  Tn;  it  is  denoted  f3n  =  /3en_1:en, 
highlighting  the  assumption  that  semi-Markov  jump  distributions  depend  only  on 
the  previous  and  current  market  states.  The  density  function  of  fyj  is  b(-  \  i,j)  and 
Let  is  defined  in  (3.2). 


Remark  4.1.  A  few  observations  about  the  model  in  the  context  of  Remark  3.2 
are  in  order.  The  solution  of  (4.1)  can  be  described  by  the  following  discrete  time 
iterative  process  Ladde  &  Ladde  (2013) 


/ 

x(t,Tn,  xn)  =  xn  exp 

/  dlfr 

< 

JTn  J 

k xn  — (5nx{Tn  ,  Tn—\ ,  xn—\ ) , 

t  G 


[Xn,  Tn_j_i) 


(4.2) 
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where  Let  is  defined  as  the  exponent  of  the  solution  process  of  (3.1)  as  expressed 
in  (3.25).  The  semi-Markov  process  decomposes  both  the  time  and  state  domains 
causing  structural  changes  in  the  stock  price  process,  while  the  Levy  process  directly 
decomposes  the  state  domain  of  definition  of  the  stochastic  dynamic  model. 

Remark  4.2.  From  (3.25)  and  (4.1),  the  size  of  the  jump  in  log  price  at  time  Tn 
is  In ((3n).  The  density  function  of  In (/3n)  is  described  by 

b(z  |  0n— l,  0n)  =  b(ez  |  0n_i,  0n)ez,  (4.3) 

where  &(•  |  0n_i,  0n)  denotes  the  density  of  /3n  and  e  is  the  Naperian  base.  We  further 
note  that  the  discrete  time  dynamic  system  in  (4.2)  is  an  intervention  process.  A 
feature  of  interest  of  this  model  is  its  potential  to  capture,  simultaneously,  three 
important  stylized  facts.  The  volatility  clustering  exhibited  in  log  return  time  series, 
the  slowly  decaying  autocorrelation  of  square  returns  and  the  observed  correlation 
between  log  returns  and  volatility  (Chourdakis  2005,  Bulla  2006).  As  the  market 
switches  from  one  state  to  another,  the  diffusion  rate  changes  while  the  asset  price 
is  subjected  to  a  jump.  Thus  the  diffusion  rate  and  the  price  jumps  are  modulated 
by  the  process  0. 

For  the  development  of  an  infinitesimal  generator,  in  the  following,  we  define 
a  point  process  encoding  both  the  regime  switches  and  the  jumps  of  x  at  regime 
switches.  At  each  regime  change,  we  note  that  the  jump  in  log  price  is  In (/3n).  We 
define  E 2  =  {(i,  j),  (i,j)  G  E  x  E,  i  ^  j}  and  the  power  set  of  E2 ,  V{E2).  B(M)  is 
the  Borel  sigma  algebra  of  the  real  line.  We  are  ready  to  define  the  aforementioned 
point  process. 

Definition  4.2.  (3n  and  6n  are  introduced  in  Definition  4.1.  Let  N(t,A,B)  be  a 
stochastic  process  defined  on  [0,T]  x  B(M)  x  P(E2)  into  (0,  oo)  as 

N(t,A,B)  =  Y:  1(t>Tri,ln(/3n)eA,(6>n_i,6>ri)eB);  (4.4) 

n>  1 

and  N(t,  A ,  B )  stands  for  the  number  of  regime  switches  in  B  with  corresponding 
log  price  jumps  ln(/?n)  G  A  by  time  t. 

Remark  4.3.  We  observe  that 

N(t,A,B)  =  ]T  (4.5) 

where  N(t,A ,  {(i,  j)})  counts  the  number  of  regime  switches  from  i  to  j  with  cor¬ 
responding  log  price  jump  In (/3n)  G  A. 

In  the  following  lemma,  we  derive  the  predictable  compensator  process  for 
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Lemma  4.1.  Let  N(t ,  A,  {( i,j )})  be  the  point  process  introduced  in  Definition  4.2. 
Then 


N(t  A  Tn,  A,  {(i,  j)})  -  7 (t  A  Tn,A ,  {(i,j)}) 
is  a  martingale  with  respect  to  the  filtration  (Ht)t>o5Vn  E  7(1,  oo),  where 


(4.6) 


7 (t,A,{(i,j)})=  f  [  b(z  |  i,j)Xij(ys)dzds,  (4.7) 

7o 


and  X ij  are  defined  in  (2.9). 


Proof.  From  Bremaud  (1981),  Siu  &  Ladde  (2011),  it  is  enough  to  prove  that 
N(t  A  Tn,  A,  {(i,  j)})  -  7 (t  A  Tn,  A,  {( i,j )})  is  an  (Hs)  ^-martingale.  For  any  0  < 
s  <  t  and  for  each  n  E  7(1,  oo),  it  satisfies 

E[[N(t  A  Tn,  A,  {(i,  j )})  -  7(t  A  Tn,  A,  {(7,  j)})] 

-  [N(s  A  Tn,  A,  {(*,  j)})  -  7(5  A  Tn,  A,  {(i,  j)})]  |  Ha]  =  0  (4.8) 

and  if  and  only  if 

E(N(t  A  Tn,  A,  {(i,  j )})  -  JV(*  A  Tn,  A,  {(*,  j )})  |  Hs) 

=  E(^(t  A  Tn,  A,  {(ij)})  -  7 (^  A  Tn,  A,  {(i,  j)})  |  fla).  (4.9) 

We  prove  that  (4.9)  holds.  We  first  prove  that  (4.9)  holds  when  the  jump  process 
N  is  stopped  at  T\.  We  then  prove  by  the  Principle  of  Mathematical  Induction 
that  (4.9)  is  true  when  N  is  stopped  at  time  Tn.  From  Definition  4.2,  (2.5)  and  for 
0  <  s  <  t,  we  have 


E(N(t  A  Ti,  A,  {(i,j)})  -  N(s  AT,,  A,  {( i,j )})  |  fls) 

^h(Ti<t,in(0i)eA,$itmj,eo=i)  -  l(Ti<s,ln(/3i)eA,0i=j,0o=i)  I  #0,Tl  >  s), 

for  T,  >  s 

0,  for  T i  <  s, 

=  l(Ti>s)^(l(T1<t,ln(/3i)eA,0i=j.0o=i)  “  l(Ti<s,ln(0i)€A,0i=.7,0o=i)  I  ^0 ,  >  s) 

=  l(Ti>s)-E'(l(s<Ti<t,ln(/3i)eA,6»i=j,6»0=j)  I  #0,-Tl  >  «) 

P(s  <  Pi  <  t,  ln(/3i)  e  A,6,  =  j\0o  =  i) 


~  l(Ti>s)l(0o=») 


-  l(T1>s)l((?o=i)Pii 


P(Ti  >  s  I  00  =  i) 

P(s  <T,M^)  G  A\e,  =  j,e0  =  i) 


S(s\i) 

P(t<T1,]n(A)eA\e1=j,0o  =  i) 
S(s\i) 
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l{Tl>s)l(90=i)P(HPi)  e  A  |  \n((3i)i,j)pij 


P{s<Ti\e0  =  i,91=j) 
S{s\i) 

P(t  <  Ti  |6>0  =  =  j) 

S'(sli) 


=  l(T1>5)l(eo=i)Pijp(ln(/3i)  e  A  |  <90  =  i,6>i  =  j)  (4-10) 

where  A S(t  \  i,j)  =  S(t\6o  =  i,0\  =  j)  —  S(s\0o  =  i,  6\  =  j),  with  5(-  |  i,  j)  denot¬ 
ing  the  conditional  survival  distribution  of  sojourn  time  when  the  process  switches 
from  i  to  j.  From  (2.8),  (2.9)  and  (4.3),  we  have 

E[N(t  A  Ti,  A,  -  N(s  A  Ti,  A,  {ij})} 

=  [  [  hT1>s)^{e0=i)Kz\i,j)Xi^j{yu)dudz. 

J  S  J  Z(zA 

On  the  other  hand,  from  (4.7)  and  (2.8),  we  obtain 
E[l{t  A  Ti,  A,  {(*,  j)})  -  7(«  A  Ti,  A,  {(*,  j)})  |  Hs] 


(4.11) 


-  l(0o=*)l n>sE 


•t  ATi 


-A  J  o 


Kz  |  i,j)Xij(yu)dudz 


sATi 


zeAJ  o 


6(2  |  i,j)\ij(yu)dudz  |  Hs 


=  lTi>sl(6»o=i)-®’ 


[TiM 

/  P(ln(/3i)  G  A  |  i,  j)Xij  (yu)du  |  Ti  >  s,  6»0  =  i 

J  Ti  As 


1^1  >S  l(0o=z) 


=  lTi>s  1(00=0 


(Fubini’s  theorem) 
-P(ln(/3i)  gA|i,i)VAt/M)dM~^f(^°_=  0 

/s  J  T\  As  b{S\Vo  —  l) 


oo  pTiAt 


nr  At 

P{ ln(/?i)  G  A  I  i,j)Xij(yu)di 

As 


—dS{r  \Qq=i) 

‘  WWo  =  i) 


rrM  P(ln(/?i)  G  A  I  i,j)Xij(yu)du  dS^r\0^ 


t  J  rAs 


l^i  >s  l(0o=O 


50 


5(spo  =  »). 

b)//v"  n(/?i)  G  A  |  i,j)Xij(yu)dudS(r  \  90) 


/OO  I"  rt 

J  P( ln(/3i)  G  A  I  i,j)Xij(yu)du 


+ 


-dS(r  |  6>0  =i) 


S(s  \0O  =i) 
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^-Ti>s^-(Oo=i) 


S(s  |  60  = 


n(/?i)  G  A\i,j)Xij(yu)dS(r\0o)du 


+ 


n 


—dS{r  |  #o  =  i) 


n  S(s  \0q  =  i) 


P{ ln(/3i)  G  A\iJ)\id(yu)du 


—  1ti>s1(0o=«)  |  0q  _  *  j  1)  £=  A  |  i,j)K,j(yu) 

f  dS(r\0o=  i)  du 

J  u 

s(t  I  00  — *)  f  p(\n^  e  A  |  ij)\itj(u)dv 
—  l)  J[s,t] 


<5(s  |  0O 

=  lT1>sl(0o=i)-P(ln(/?i)  G  4  |  i,j) 

S(t\0o  =  i)  rt 


•  S(i  |  0O  =  *)]  du  +  /  Vi  (Vu)du 

O  (5  |  C7o  —  t)  J s 


-  lTi>sl(0o=i) 


S(s  |  0O  = 
P(ln(^)ei4|i,i) 


=  l-T1>sl(80=i)Pij 


S(S  |  0Q  =  *) 

P(ln(/3i)  G  ^4  |  i,  j)  rt 
S(a  |  0o  =  i) 


/  Xij(u)S(u  |  #o  = 

J  s 

ha,  |  n  dS(u\i,j ) 

J  S{u  1 0o  =  *)- 


S(u  |  0o  =  *) 


(4.12) 


=  Iti  >s  1  (6io=i)4>ii-f>0n( A  )  £  V 
From  (4.10)  and  (4.12),  we  get 

£(7(*  A  Ti,  A  {(*,  j)})  -  7 (s  A  Ti,  A  {(*,  j)})  |  fls) 

=  w/jo  g  a  i 

=  E(N(t  ATi,A,  {( i,j )})  -  IV  (s  A  T1;  A,  {(i,j)})  |  Hs). 

This  establishes  (4.9).  Hence,  the  stopped  point  process  N(t  A  Ti,  A,  has 

predictable  compensator  j(t  ATi,  A  x  {(i,j)})  defined  in  (4.7).  Assuming  that  (4.8) 
is  valid  for  some  k  G  7(1,  oo),  and  repeating  the  above  argument,  we  verify  the 
induction  assumption.  By  the  principle  of  mathematical  induction,  we  conclude 
that  N((t  ATk,t  ATk+i],A,{(i,j)})  -  A  Tk,t  A  Tk+1],  A,  {(i,  j)})  is  an  (Ht)t>o- 
martingale.  □ 


Prior  to  turning  our  attention  to  the  infinitesimal  generator,  we  first  establish 
Ito  differential  formula  for  (4.1). 
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Theorem  4.1  (Ito  Differential  Formula).  Let  V  G  C[M+  x  M+  x  M+  x  M,  M]  be 

continuously  differentiable  in  the  first  and  second  variables  and  twice  continuously 
differentiable  function  in  the  fourth  variable.  Let  x ,  y ,  N  and  7  be  stochastic  pro¬ 
cesses  defined  in  (4.1),  (2.14),  (4.4)  and  (4.7),  respectively.  Moreover ,  processes  N 
and  if  do  not  jump  simultaneously  P -almost  surely.  Then 

dV(s,ys-,Qs-,xs-) 

dV 

(£F)(s,  ys-  1  ds- ,  xs—  )ds  H-  cr{6s—  )xs—  dBs 


+ 


+ 


/  [V(s,ys,Qs,xs-  +xs-G(z,Qs))  -  V(ys,6s,xs-)\fi(6s,dz,ds) 
J\z\<l 

/  [V(s,ys,Qg,xs-  +Xs-H(z,6s ))  -  V(s,  ys,  Qs,  xs-)]ip(0s:  dz,  ds) 

J\z\>l 

+  Y  {V(s>ysJ,xs-ez) 

JzeR  7t> \  r/)  i 

(4.13) 

d2V 


k i€E\{»s -} 

-  ^ (s,  2/s- ,  0s- ,  Zs-  )}iV'(ds,  {(<9S- ,  /)}), 

/or  0a-  E  E,  where 

dV  1  2//i  n  2 


dV 

dV 

ds 

+  ( 

4^(5,  ys~  •> 

J\z\<l 

-V{s. 

)  2/s  5  0S  >  ^S 

+  l 

^(S,  J/s-,( 

J\z\>l 

-V(s, 

y8 

0  T 

-  ,  Vg-  ,  Ujg- 

7, 

M 

1 

'3&s- 

ay 


v(Os-,dz) 


'  V (s,  ys-  j  6S—  j  xs- )]  b(z  |  $s—  ,  j)dz. 


(4.14) 


0S-  G  E  and  N  =  TV  —  7. 


Proof.  Let  V  be  defined  as  in  the  theorem.  Let  {Tn}ff=1  be  a  sequence  of  semi- 
Markov  jump  times  and  To  =  0.  For  t  G  M+,  we  can  find  an  interval  [Tn,Tn+i] 
such  that  Tn  <  t  <  t  +  At  <  Tn+ 1  for  some  n  G  N.  Let  {Jj1}*?” 0  C  [Tn,Tn+i] 
and  J0n  =  Tn  be  a  finite  sequence  of  jump  times  due  to  the  Levy  jump  process  for 
kn  G  N.  We  further  note  that  the  interval  can  be  rewritten  as 

[Tn,  Tn+ 1]  =  [Tn,  T~+1]  U  [T~+1,  Tn+l}.  (4.15) 
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We  observe  that  [Jj\  J™+1\  fl  [JJY 1,  Jj'+i]  =  0.  In  addition, 

kn 

[Tn,  T-+ 1]  =  JjVl]  u  v;+i]).  (4.16) 

3= 0 

It  is  known  that  the  state  dynamic  process  operating  under  the  above  stated 
conditions  decomposes  into  three  parts,  namely,  the  continuous  time,  the  Levy  jump 
time  and  the  semi-Markov  jump  time.  In  fact,  the  solution  process  of  (4.1) / (3.1) 
can  be  rewritten  as 


xt  =  xct  +  xd  +  xst, 


(4.17) 


where  xd  and  xf  are  due  to  the  presence  of  continuous  process,  Levy  process 
and  semi-Markov  process,  respectively.  We  further  observe  that  for  s  G  [Tn,Tn+ 1], 
we  have:  s  =  s~  +  (s  —  s~)  =  s~  +  As,  where  As  =  s  —  s~ ,  s~  ^  s.  From 
Definitions  2.3,  4.1,  we  note  that  ys  =  ys-  and  0S  =  0S-  for  s  G  [Tn,  T-+ 1]  and  for 
5  =  Tn+ 1,  s  ^  s~,  ys-  ±  ys  and  <9Tn+1  /  0S- .  Moreover,  there  is  a  j  G  7(1,  kn  -  1) 
such  that  s  G  [Jj\  J  U  [J%n,T~+1\.  We  choose  As  so  that  s  +  As  G  [Jj\  J™+1]. 
For  these  choices  of  s  and  s  +  As,  we  have 


Z/s+As  —  Vs~  “1“ 

<  =  @ s~ 

Xs+As  =xs-  +  Ax*, 

Furthermore, 


'  Axca 


Axds , 
Ax*, 


ifsG[j;,JJVi)U[J-n,Tn+i), 

for  j  G  7(0,  kn  —  1)  and  n  G  7(0,  oc) 
if  s  =  JJVi,  for  j  G  7(1,  fcn  -  1) 
if  s  =  Tn+ 1 


(4.18) 


(4.19) 


(4.19)  implies  that  for  s,  As  G  [Jj\  JJYi),  the  change  in  state  dynamic  process  is  in 
the  absence  of  the  influence  of  Levy  jump  process.  On  the  other  hand,  for  s  =  Jj1  for 
each  j  G  7(1,  kn),  the  dynamic  process  is  interrupted  by  the  presence  of  Levy  jumps. 
Finally,  if  s  =  Tn+ j,  then  the  dynamic  system  undergoes  a  structural  change.  Here 
the  structural  change  is  under  the  influence  of  the  semi-Markov  process.  Therefore, 
there  is  no  contribution  of  the  continuous  time  dynamic  process. 

Based  on  the  nature  of  the  dynamic  process  operating  under  continuous  time 
process,  semi-Markov  process  and  Levy  process,  we  compute  the  change  in  the 
auxiliary  function  V  as 


V (s  H-  As,  ys+Asi  ^s+As)  I7(s,  j/s,  xs)  (4.20) 

in  the  context  of  state  dynamic  model  (3.1). 

The  computation  of  change  in  (4.20)  depends  on  the  computation  of  changes  over 
the  time  domain  of  decomposition  of  [Tn,Tn+i]  for  n  G  7(0,  oc).  For  computation 
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on  U-=o  V”,  Jj+i)  U  K  ,Tn+i],  we  utilize  the  generalized  mean  value  theorem. 
For  this  purpose,  we  pick  s,  s  +  As  G  [Jj1  ,  Jj+i)  C  [Tn,T“].  From  (4.18),  the 
computation  of  (4.20)  on  the  time  domain  Uj=o[^7’  ^j+i]  U  [^fcn ,  ^n+i]  regarding 
the  continuous  part  of  state  dynamic  (4.1) / (3.1)  is  as  follows:  The  decomposition 
of  three  subsets  of  time  domain  [Tn  ,  Tn+1],  namely,  L feP",  P+i)  U  [^B,  T~+1],  or 
JJVi],  or  [T“+,Tn+i]  for  n,  kn  G  1(0,  oo). 

V (s  “t”  As,  |/s-|-As?  ^s+As?  «^s+As)  V (^?  I/s?  ^s?  ^s) 

r1  ray 

=  J  -^-(s  +  riAs,ys+riAs,0s,xs+T]Axs)As 
dV , 


+  -{,  +  1A,,y,+yAy.,e„x.  +  yA,.)Ay, 
dV 

+  ~dx  (S  +  T]As’  Vs  +  ®s' Xs  +  VAxs)^xs 


dr] 


dV  dV  dV 

xs^)As  4  ^^-(5,  t/s,  0S,  xs^)Ays  +  (^  Z/s?  3?s)Aa?s 


l 


dx 

+  £s,y(As), 

where 

£s,y{^s) 


1  rdV  dV 

(5  +  77AS,  2/a  +  6>s,  +  ?]Ars)  -  — (s,  y8, 6S,  xs) 


dx 


Axsdr] 

(4.21) 


f1  ray  ay 

J  -g^(s  +  yAs,ys  +  r]Ays,9s,xs  +  yAxs)  -  —  (s,ys,9s,xs) 


Asdy 


+ 


f 


dV 


ay. 


(s  +  yAs,  ys  +  yA y8, 9a,xs  +  yAxa)  -  -^-(s,  ys,9s,xa) 
oy  oy 


A  ysdy. 
(4.22) 


We  again  apply  the  generalized  mean  value  theorem  to  the  integrand  in  (4.21),  and 
we  obtain 

dV  dV 

—  (s  +  yAs,ya  +  r)A ys,9s,xs  +yAxs )  -  —  (s  +  yAs,ya  +  yAya,9a,xa) 


a2y 

dx 2 

a2y 


(s,ys,9s,xs)yAxs  + 

Jo 


a2y 

ax2 


ax2 


(• s,ys,9s,xs ) 


Axsrjde. 


(s  +  As,  ys  +  rjAys ,9s,xs  +  eyAxs ) 


(4.23) 
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From  (4.23),  the  fourth  term  in  (4.21)  reduces  to 

iOT/  1 

A  xsdrj 


~0V ,  .  , 

<9w  , 

L 

—  (s  +  r?As,  ys  +  r?Ays,  03,  xs  +  r?Aa;s) 

-  -7^(s,ys,0s,xs) 

1  d2V 

=  2~dx2^S’  Vs^  @sixs)(Axs)  “I-  sx{As')^ 


(4.24) 


where 


l  ri  r 


b(A s)  =  f  f 

Jo  Jo 


d2V  d2V 

-^j(s  +  As,ys+r/Ays,es,xs  +  er/Axs)  -  -^(s,ys,9s,xs) 


+ 


f 


dV 

(s  +  rjAs,ys  +  rjAys,  6S,  xs)  -  —  (s,ys,6s,xs) 


x  ( Axs)2r]dr]de 

rl\dV_ 
dx 

From  (4.21)  and  (4.24),  we  have 

V(s  +  As,ys  +  Ays,9s  +  As  5  )  -  V(s,ys,9s,xs) 


dr]Axs . 


dV .  .  ,A  dV .  .  ,A  dV  . 

=  —  (s,ys,9s,xs)As  +  —  {s,ys,0s,xs)Ays  +  —  (s,ys,0s,xs)  Ax., 


dy 

1  d2V 

+  ^-Q^{s,ysies,xs){Axs)2  +e{As), 


(4.25) 


where  e(As)  =  sS:y(As)  +  ex(As).  The  expressions  in  (4.21),  (4.23)  and  (4.24)  are 
valid  for  all  s  G  [Jj1,  U  [  J^n ,  Tn+1]  and  for  all  j  G  7(1,  &n  —  1)  and  hence  they 
are  valid  on  the  intervals  [Tn,  Tn+i)  for  n,  kn  G  7(0,  oc). 

Using  Levy  integrals  and  a  single  jump  value,  we  compute  (4.20)  under  the 
influence  of  Levy  jump  process.  For  this  case,  we  first  compute  U(J™+1)  —  V ( Jj!^), 
where  V{J?+1)  =  V(s  A  Jj+i,y3AJ?+1,08AJ?+1,xsAj?-i  +  xsAJ^G(z,  <Wt>+1))  and 

We  set  and  compute: 

=  [V(s  A  J”+1 ,  yiA jy+1 ,  jy+1 ,  *SA Jj*+-  +  xsA j»+-  G(z,  fliA jy+1 )) 

-  r (s,  ys- ,  0S- ,  za-  )]V>(0«- ,  Az,  As) 

+  [V(s  A  J"+1 ,  Vs/\JV-+1 ,  ^.sA.7”+,  ,  ^sA J"+1 )) 

-  r (s,  y8- ,  0S- ,  za-  )]ip(03- ,  Az,  As).  (4.26) 

From  (4.26),  for  any  s  G  [Jj+i,  JJYi],  J  £  7(0,  kn  —  1)  and  n,  kn  G  7(0,  oc)  we  have 

VisAJ^-VisAJ^) 

/s+ As  /* 

/  [U(5,7/s,6>s,xs-  +XS-G(^,6>S))  -  U(5,7/s,6>s,xs-)]'0(6>s,d2;,d8) 

7m<i 
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ps-\-/S.s  p 

/  /  [V(s,ys,6s,xs-  +xs-H(z,9s))  -V(s,ys,6s,xs-)\ 

Js  J\z  |>1 


+ 


l\z\>  1 
x  ip(Os,  dz ,  ds). 


(4.27) 


The  expression  in  (4.27)  is  over  a  subinterval  U j=o  1  t^j+i ’  Jj+ l 1  °f  Pnj^n+i]- 
Finally,  for  5  G  [T“+1,Tn+i],  and  imitating  the  above  argument,  we  compute  (4.20) 
under  the  presence  of  semi-Markov  jump  as  follows: 

V(sA  -^n+ 1 5  VsATn+i  •>  @ s/\Tn-\-\  ’  XsATy+1  +  ^ XsAT~+1  ) 

-  V(s  A  Tn+1,ysAT-^,  0SAT-+11  XsAT~+1) 

=  V(s  A  Tn_|_i,  ysAT- ,  +As,  ®SAT~  .+As’  XsAT“,  n  e^)  —  2/s_  5  0 s~  3  ) 

n+1  n  +  1  n  +  1 

/•s+As 

=  /  /  [V(w,yu-,0u-,zu-ez) 

ds  d  2GM 

(4.28) 


-  r (w,  yu- ,  ,  xu-  )]iV(ds,  dz,  {0T-+i ,  #r„+1 }), 


hence,  adding  and  subtracting 

fS+As 

J  z  CM 


pS+l\S  p 

I  I  [V('U,  yu- ,  0U—  •)  xu—  c  )  V(iL)  yu— ,  0 u— ,  xu—  )]y (ds,  dz,  {0^-  5  @u\ ) 5 
ds  d^eM 

we  obtain 

V(s  A  Tn+ 1 ,  |/sATn+i  5  $sATn+i  5  XsATy+1  ^XsATy+1 ) 

—  V(s  A  Tn+1 ,  2/sAT-+i  ,  #sAT-+i  ,  xs/kT -+i  ) 

/*s+As  p 

=  /  [r(ti,2/u-,0„-,a;u-ei:)  ^y(«,yu-,0u-,a;u-)]7(du,d2;,  {0U-,0U}) 

ds  d^GM 

/>s+As  /•  ^ 

I  I  yu—  5  5  ^u_  ^  )  ^(^A  du~  5  5  )]-^V"(d,u,  dz,  ,  $u}) • 

ds  d^GM 

(4.29) 


+ 


This  expression  is  on  [Tn+1,Tn+ 1]  for  n  G  7(0,  00).  From  (4.25),  (4.27)  and  (4.29), 
(4.20)  reduces  to 


V^s  +  As^+asA+As^s+As)  -  V(s,ys,0s,xs) 

dV  dV  dV 

—  xs)A5  4  ^^-(s,  £/s,  $s,  xs)Ays  (s,  t/s,  $s,  ^s)Ars 


2  dx2 


(s,ys,^,rs)(Ars 
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rS  +  As  p 

/  /  \V Vui  tu-  H-  xu-  G(z ,  9U- )) 

Js  J\z\<  l 


+ 


'N<1 

V(u^  yu—  ,  5  )]  *0(^14)  ^^5  did) 

pS~\~As 


+ 


/>S  +  ZAS  p 

I  I  l/uj  ^U)  Ou—  )) 

Js  J|z|>l 


-  ^ (w,  2/«- ,  #u- ,  ar«- )] V’(^U)  du) 

rS+As 


+ 


/S  +  ZAS  /> 

I  ^  ^  [^(^5  ) 

JzeR  oueE\{eu -} 


^(^b  Vu~  5  ^U-  5  ^W_  )]t(^®5  djZ,  \0U-  ,  ) 

pS-\-As  p 

I  /  [V^IX,  Vu~  5  ^U_  5  ^U—  ^  ) 

J s  J z£M 

-  V (u,  yu- ,  Ou- ,  xw-  )]iV(ds,  dz,  {0U- ,  0U})  +  e(A s). 


(4.30) 


For  small  As,  applying  the  concepts  of  stochastic  differentials  Applebaum  (2009), 
adding  and  subtracting 


rs+As  p 

I  I  \Y  (^?  Vu-)  @u-> 

Js  J\z\>l 

V {u,  yu- ,  0U- ,  xu- 


/»S  +  As  p 

I  I  \y  (^?  Hu')  @W) 

Js  J\z\<l 

V(lLl  yu—  ,  ^n_  5  Ty  — 


Xu-  +xu-H(z,e  „-)) 
)]z/(0„- ,  dz)du  and 
-|-  Xn—  Cr(z,  )) 
)]z/(6>M-,dz)dw, 


(4.30)  reduces  to 

dV 

dV(s,ys-,9s-,xs~)  =  CV(s,ys-,9s-,xs-)ds  +  a(9s-)xs-  —(s,ys,9s,xs)dBs 

+  /  [V(s,ys,9s,xs-  +xs-H(z,9s-)xs-) 

J\z\>l 

-V(s,ys-,9s-,xs-)]'ip(9s,dz,ds) 

+  /  [V(s,ys,9s,xs-  +  xs-G(z,9s-)xs~) 

J\z\<l 

-  V (s,  ys- ,  9S- ,  xs-  )}ip{9s,  dz,  ds) 

+  /  [ V(s,ys-,9s-,xs-ez ) 

JzeR 

~V(s,ys-,9s-,xs-)]N(ds,dz  x  {0S-,0S}).  (4.31) 
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This  establishes  Ito  differential  formula  (4.13)  for  Levy  type  stochastic  differential 
equation  under  semi- Markovian  structural  perturbations.  Here  C  in  (4.14)  is  the 
linear  differential  operator  relative  to  (4.1).  □ 


In  the  following,  based  on  Theorem  4.1,  we  present  a  concept  of  infinitesimal 
generator  and  a  few  results  as  special  cases. 


Definition  4.3.  For  the  function  V  defined  in  Theorem  4.1  and  using  (4.30),  an 
infinitesimal  generator  of  (4.1)  is  defined  by 


lim 

At— s-0 


At 


E[V (t  T  At,  2/t+Atj  ^t+At?  *^t+At) 


-V(t,yt,etixt)  |  yt  =  Vt-,0t  =  0t-,xt  =  xt~] 

=  AV (t,yt-,0t-,xt-),  for  0t-  e  E.  (4.32) 

Moreover,  a  one  parameter  family  of  semi-group  is  generated  by 
dV 

-Qjr(t,yt-,6t-,  xt- )  =  AV(t,yt-,0t-,  xt- ) ,  (4.33) 

where  A  =  C  in  (4.13)  and  ^r(t,yt,0t,xt)  is  the  conditional  partial  derivative 
defined  by  the  left-hand  side  of  (4.32). 


We  present  special  cases  of  the  developed  infinitesimal  generator  in  Defini¬ 
tion  4.3. 


Remark  4.4.  From  Remark  3.2,  the  infinitesimal  generator  A  defined  in  (4.32) 
extends  the  earlier  work  in  a  systematic  way.  In  fact,  this  generator  includes 
the  infinitesimal  generator  influenced  by  finite  state  Markov  chain  (Elliott  et  al 
2005,  Chourdakis  2005).  Moreover,  it  also  includes  the  generator  influenced  by  a 
finite  state  semi-Markov  process  (Ghosh  &  Goswami  2009,  Hunt  &  Hahn  2010, 
Swishchuk  &  Islam  2011).  If  H  and  G  are  replaced  by  eG  —  1  and  eH  —  1,  then  A 
in  (4.32)  in  the  context  of  (4.14)  is 


v4E(«s,  ys— ,  0  s— ,  0S ,  xs— ) 

3V  3V  ✓  „  v  dV  1  2  2sn  ^d2V 
-^  +  ^  +  Kds-)Xs-d^  +  2Xs~(j{es-)dP 


+ 


/ 

J\z\<l 


V(s,ys-,6s-,6s,xs-  +xs-[eG<'z’e*-^  -  1]) 

dV 


V{s,ys-,Qs-,Qs,xs-)  —  xs-  [eG^,e*-)  -  1] 


dx 


v{0s- ,  dz) 


+  f  [V («,  Vs-  ■,  0S- ,  0s,xs-  +  xs-  }  -  1]) 

J\z\>l 
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^(^5  Vs~  i  d  s~  5  ®si  Xs—  i  dz^j 

+  ^2  ^es-gs(ysW(s,ys-,Os,xs-ez) 

JzER  eseE,es^Gs- 

~  V (s,  Vs- ,  Os- ,  xs-  )\KZ  I  Os- ,  Os)dz. 

A  few  notes  on  the  nature  of  the  infinitesimal  operator. 

Remark  4.5.  We  further  remark  that  the  infinitesimal  generator  defined  in  (4.32) 
can  be  rewritten  in  a  m  x  m  matrix  form.  In  fact  the  partial  differential  equa¬ 
tions  in  (4.33)  are  a  system  of  partial  differential  equations  of  dimension  m.  More 
precisely,  (4.33)  is  a  linear  system  of  partial  differential  equations  with  variable 
coefficients. 


Remark  4.6.  For  V(t,ytl0t,xt)  =  xt,  the  conclusion  of  Theorem  4.1  reduces  to 

dxt  =  CV(t,yt-,6t-,xt-)dt  +  a(d^)dBt  +  f  xt-G(z,et-)ip(6t-,dz,dt) 

J\z\<l 

+  f  xt-H(z,0t-)$(0t-,dz,dt) 

J\z\>l 

+  [  [xt-ieZ -l)N{dt,dz,{{Ot-,j)})\,  (4.34) 

JzeR  jeE\{6t-} 

where 


with 


cv  (t,  vt- ,  ot~ ,  xt~ ) 


xt- 


y{0t- )  +  f  H(z,  0t- )v{0t- ,  dz) 
J\z\>l 


+ 


I  [(ez  -  l)*et-,j(yt-)]Kz\Ot-,j)dz 

J\z\>i  jeE\{et_} 


dxt  =  xt-dM ): 


(4.35) 

(4.36) 


dM°  = 


y(0t-  )dt  +  f  G(z,Ot-)m- 
+  [  H(z,6t-)ip(6t-,dz,dt) 


,  dz ,  dt) 


J\z\>l 

+  [  H  [(ez-l  )N(dt,dz,{(8t-,j)})] 

J  zeR  jeE\{9t-} 
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+  Y  l(eZ -1)Xot-d(yt-)]b(z\Ot-,j)dzdt 

jeE\{et-} 


+  /  H(Ot- ,  z)v(Qt-  ,dz)dt 
J\z\>l 

=  dLet+  [  Y  [(e*-l  )N(dt,dz,{(Ot-,j)})\,  (4.37) 

JzeRjeE\{9t _} 

where  Le  is  defined  in  (3.2).  Furthermore,  we  note  that  the  solution  process  deter- 
mined  by  (4.36)  has  another  solution  representation  of  (4.1)  in  the  framework  of 
Remark  4.2.  In  fact,  the  closed  form  solution  representation  of  (4.36)  is  as  follows: 


xt  =  xq  exp 


+ 


[  [  J2  [ zN{ds,dz,{(6s-,j )})] 

Jo  Jze ®jeE\{es_} 

f  n(0s-)ds-\  [  a2(ds-)ds 
Jo  z  Jo 


+  f  f  [ln(l  +  G(z,6s-))  -  G(z,9s-)]u(9s-,dz)ds 
Jo  J|«|<i 

+  [  cr{Os-)dBs+  [  [  ln(l  +  G(z,0s-))tp(0s- ,dz,ds) 
Jo  Jo  j\z\<l 

+  f  f  ln(l  +  H(z,  6S-  ))^{0S- ,  dz,  ds)  . 

JO  J\z\>l 


(4.38) 


In  the  following  section,  we  utilize  the  infinitesimal  generator  of  the  exponen¬ 
tial  semi-Markov  Levy  switching  process  to  find  a  closed  form  expression  of  the 
characteristic  function  of  the  ln(xt)  from  (4.1). 


5.  Characteristic  Function 


In  this  section,  we  derive  a  closed  form  expression  for  the  conditional  characteristic 
function 


■*(u,t,y,j,x)=E[eiul^\y0 


y,  80  =  j,  x0  =  x } 


(5.1) 


of  the  log  price  process 


n(t) 

\n(xt)  =  YHPp)  +  Let , 

P=  1 


(5.2) 


where  /?n,  L°  in  Definition  4.1  and  xt  is  the  closed  form  solution  process  of  (4.1)  in 
the  context  of  (3.25). 


Lemma  5.1.  Let  Let,  x ,  y  and  7  be  defined  in  (3.2),  (4.1),  (2.14)  and  (4.7),  respec¬ 
tively.  A  closed  form  expression  for  the  conditional  characteristic  vector  function  of 
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ln(x)  is 


^!(u,t,y,x)  =  ex.p[iu  ln(x)]  exp 


'  rt+y 

/  M(u,  s)ds 

Jy 


1, 


(5.3) 


where  i  =  \/ — 1;  £,  y ,  x)  zs  an  m-dimensional  column  vector  with  kth  component 

4/(iq  £,  y,  &,  x),  for  k  £  E;  1  is  m  x  1  vector  with  components  ones ,  and  M(u,y)  = 
(MPg(w,t/))mxm  is  an  m  x  m  matrix  defined  by 


y )  —  < 


iuy(q)  —  -a2(q)u2  +  /  —  1  —  iuG(z,  q)]v(q,  dz) 

2  ./|*|<i 

+  [  [eiuH^  -  1]%,  dz)  +  \q:q(y),  if  p  =  q 

J\z  |>1 


A q,p(y)  /  eluzb(z\q,p)dz ,  otherwise , 

and  it  is  assumed  to  satisfy  the  Lie  bracket-type  condition 

[M(u,y1),M(u,y2)]=0,  Vyi,y2eM+. 


(5.4) 


Proof.  From  (5.3),  first  we  observe  that 

\k(iq  £,  z/,  x)  =  exp[mL^]  exp 


r 

/  M(rq  s)d6 

V  ?v 


1. 


(5.5) 


We  note  that  T(iq  £,  y,  0tl  x)  possesses  all  smoothness  properties  of  V  defined  in 
Theorem  4.1.  Therefore,  following  the  argument  used  in  the  proof  of  Theorem  4.1, 
Definition  4.3  and  Remark  4.5,  we  conclude  that  T  is  in  the  domain  of  the  infinites¬ 
imal  generator  of  the  process  (yt,  0t,  iuL^)te^  Ty  Moreover,  it  satisfies  the  following 
system  of  linear  partial  differential  equation: 


&$(u,t,y,k,  x)  =A t  k  )  hvkeE,  (5.6) 

Ot 

where  A  is  the  operator  defined  in  Definition  4.3.  From  Remark  4.4  with  y(0s-) 
replaced  by 

A»(0a- )  +  \X0e- )  +  f  [eGiz’e°- }  -  G(Z,  0S- )  -  1  ]i/(08- ,  dz) 

Z  J\z\<l 


and  for  T  defined  in  (5.1),  we  have 


AAt (u^  5,  ys— ,  0 s— ,  x) 

d^f(u,  s,  y,  6S- ,  x)  d^(u,s,y,6s-,x) 


ds 


m(^)  +  2cr2(^) 


dy 

(pHup^,  y,9s-,x) 

dx 
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+ 


+ 


/ 

^(w,s 

>  y  1 9  s 

-,x  +  x[eiuG{- 

J\z\<l 

'Uw,  S. 

i  V : 9 s~ ' 

,x)  - 

iuxG(z ,  0S-)- 

/  [*(«,«, 

y,9s- 

)X  +  a;[ei«^ 

•'M>1 

°-i]) 

d^(u,s,y,6s-,x) 


dx 


v(6s- ,  dz) 


i  q  \i  (r.  j  \  1  2  2  2 in  \ ^  Ui  9  s  —  ?  *^) 

-^(u,s,y,9s-,x)\v(9s-,dz)  -  -xsu  a  (9S~) - - 


+ 


f  J2  *(%s,y,Os,xeiuz)b(z\es-,es)\os_,es(y)dz 
JzeR  e3e E\{es-} 


-'&(u,s,y,08-,x)\e  g  (y)  for  9S-  e  E. 


(5.7) 


Now,  we  assume  that  ^(u,  £,  y,  &,  x)  =  exp[i?x  ln(x)]h(?i,  £,  y,  fc),  where  h(u,t,y,k) 
is  the  /cth  component  of  an  unknown  m- dimensional  vector  function  h(iq  £,  y)  = 
[ft(tq  £,  y,  1), . . . ,  h(u,  t ,  t/,  m)]T.  From  this,  (5.6)  reduces  to  the  following  system  of 
partial  differential  equations: 


dh(u,t,y,k)  dh(u,t,y,k) 


dt 


dy 


+  h(u,t,y,k) 


y(k)  +  2<j2(^) 


+  \o-2(k)[—iu  —  u2]  +  f  ^uG(z,k)  _  i  _  iuG(z,  k)]u(k,  dz) 
2  7|2|<i 


+  f  [eiuH^  -l]v(k,dz)  +  \ktk(y) 

J\z\>l 

+  f  'P\k,j(y)h(u,t,y,j)(y)eluzb(z\k,j)dz 
JzeR^k 


dh(u ,  t,  y,  k) 
dy 


+  h(u,t,y,k) 


iufjb(k) 


\k)u 2 


+  /  [eiwGf(*’fe)  -  1  -  iuG(z,  fc)]i/(fc,  dz) 
dp|<i 

+  [  [eiuH^  -lMk,dz)  +  \ktk(y) 

J\z\>l 

+  [  J2  xk,j  (y)h(u,  t,  y,  j) (y)eluzb(z  \  k,  j)dz. 


(5.8) 


As  stated  in  Remark  4.5,  the  coefficients  of  h  are  defined  by  the  elements  of  A  asso¬ 
ciated  with  £,  y ,  x),  in  particular,  the  mxm  matrix  M(u,  y)  =  y))mxm 
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defined  in  (5.4).  From  the  definition  of  h (u,t,y),  (5.8)  reduces  to 
dh(u,t,y)  _  dh(u,t,y) 


d  t 


dy 


+  M(u,y)h(u,t,y),  h(u,  0,  y)  =  1  =  (1, . . . ,  1)  . 


(5.9) 


Using  the  method  of  characteristics,  the  system  of  partial  differential  equations  (5.9) 

dy 
dt 


can  be  solved.  In  this  case,  the  characteristic  curves  are  determined  by  =  ±1. 


Solving  these  differential  equations,  we  obtain 

7]  =  t  —  y  and  £  =  t  +  y.  (5.10) 

We  use  the  above  change  of  variable  to  define  the  transforms  h  and  M  from  h  and 
M,  respectively,  as  functions  of  (rjX) 

77 +  C  -V  +  C 


=  hU, 


C)  =M[u, 


2  ’  2 
-77  +  C 


(5.11) 


From  (5.11),  the  initial  value  problem  (5.9)  reduces  to  the  ODE 

dh^dr]'  ^  =  ^  ^  Hu,  -y, V )  =  1  (5-12) 

Under  condition  (5.4),  the  general  solution  of  the  linear  homogeneous  ODE  with 
time  varying  coefficients  is  Magnus  (1954) 

rp 


h(7 7,77,0  =  exP 


1  n  ~ 

2  J  M(u,K,()d* 


•9(0, 


(5.13) 


where  g  is  an  arbitrary  m- dimensional  vector  function.  Using  the  initial  condition 
in  (5.12),  g  is  determined  by 

ro 


g{  0  =  exp 


1  r  ~ 

2  J 


VCe  [o,T], 


This  together  with  (5.13),  yields  the  solution  of  the  initial  value  problem  (5.12)  as 


h(u,77,C)  =  exP 


1  n  „ 

2  J  M(u,K,C)dK 


Using  the  inverse  of  the  transformation  defined  in  (5.10),  the  solution  of  the  original 
initial  value  problem  (5.9)  becomes 

ft-y 

h (u,  t,  y )  =  exp 

'  '-t-y 


-  exp 


1 

2  / 

Z  d-t-y 
"  rt+y 

/  M(u ,  s)ds 

-dy 


-K+t+y  ' 
u, - - - )cLk, 


This  establishes  the  conditional  characteristic  function  for  the  log  prices  described 
by  (4.1).  □ 
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Remark  5.1.  We  note  that  the  closed  form  exponential  expression  (5.3)  holds  only 
under  condition  (5.4).  This  is  due  to  the  fact  that  system  of  ode  in  (5.12)  has  time 
varying  coefficients.  Assuming  the  matrix  (Mij)mxrn  has  continuous  entries  with 
respect  to  y ,  the  system  of  differential  equations  with  time  varying  coefficients  (5.12) 
has  a  fundamental  matrix  <£(?/,  t,  i/).  Therefore,  the  characteristic  function  in  (5.1) 
is  described  by 

At(u,t,y,x)  =  exp[m  ln(x)]  [<L(?i,  £,  y)Q(u,  0,  y)_1]l.  (5-14) 

The  characteristic  function  T(rq  £,  ?/,  x)  has  a  closed  form  expression  if  4>  has  a 
closed  form  expression.  Theorem  5.1  corresponds  to  the  particular  case  where 
V(u,t,y,x)  =  exp [f*+v  M(u,s)ds}. 

In  the  following,  we  present  a  consequence  of  Lemma  5.1  that  extends  the  char¬ 
acteristic  function  of  sojourn  time  of  finite  state  Markov  processes  Momeya  (2012) 
and  Elliott  &  Osakwe  (2006). 


Corollary  5.1.  We  denote  01(h)  the  time  spent  by  the  semi-Markov  process 
(0t)te[o,T]  in  its  state  k  in  the  time  interval  [s, t],  Vs,f  G  [0,T],  with  s  <  £,  and 
0\  =  \Ol(  1),  0^(2), ... ,  0\(m)]T  denotes  the  m-dimensional  occupation  time  vector 
of  6.  If  a  =  (ai,  a2, . .  • ,  am)T  is  an  m  x  1  vector  of  constant  real  numbers ,  then 


where 


E[eiu^\ys=y,es=k\ 


exp 


M(u ,  s)ds 


Mp,q(u,y) 


iuaq  +  Xqtq(y),  if  p  =  q 
A  p,q(y),  otherwise. 


(5.15) 


Proof.  Let  (dt)te[o,T]  denote  a  stochastic  process  with  at  =  aj  whenever  0n^  =  j. 
exp [m(0*,a)]  =  exp  iu  akOl(k )  ,  for  i  =  \/— 1 


=  exp 


=  exp 


iu  E  akOl{k) 

.  keE 
't 

iu  I  a 


.  f 

IU  / 

J  s 

.  [ 

m  / 

J  s 


-dv 


,  (as  a  is  piecewise  constant) 


iu  /  av-dv 


=  elufsdLv^  with  Lev  =  av,  Vr  E  [0,T]. 

The  characteristic  function  of  the  semi-Markov  occupation  times  in  the  time  interval 
[s,t\  becomes 

E[eiu^\ys=y,9s  =  k] 

=  E[eiuVdL"\ys=y,es  =  k } 

=  E[eiuL°-s  1 2/0  =  j/,  0O  =  &], 
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since  the  couple  (0,y)  is  homogeneous.  Applying  Lemma  5.1  with  /3n  =  1  Vn  G 
7(0,  oo)  proves  the  result.  Corollary  5.1  is  a  direct  extension  of  the  results  in  Hain- 
aut  &  Colwell  (2014),  Buffington  &  Elliott  (2002)  and  Momeya  (2012).  □ 


6.  Change  of  Measure  and  Pricing  Kernels 

In  this  section,  we  introduce  the  conditional  minimum  equivalent  martingale  mea¬ 
suring  the  semi-Markov  jump  risk  and  the  Levy  risk.  In  addition,  we  also  develop 
an  unconditional  minimum  entropy  martingale  measure  and  the  Esscher  transform 
Siu  &  Yang  (2009)  measuring  all  three  risks,  namely,  Levy  risk,  semi-Markov  jump 
risk  and  regime  switching  risk.  Prior  to  the  development  of  these  concepts  and 
results,  we  utilize  the  closed  form  solution  representation  of  (3.1)  to  shed  a  light  on 
the  martingale  property  of  the  solution  process  of  the  Levy  type  stochastic  linear 
differential  equations.  For  this  purpose,  let  xt  be  the  solution  process  of  (3.1)  and 
assume  that  it  is  a  (Mt  V  Lt)- mart  ingale,  that  is,  for  s  <  £,  E[xt  —  xs  |  Hs  V  Ls]  =  0. 
This  is  represented  by  the  following  illustrations. 

Illustration  6.1. 


(1)  From  E[xt  —  xs  |  Hs  V  Ls]  =  0  and  (3.25),  it  is  obvious  that  the  solution  process 
xt  of  (3.1) 

pt  pt  pt 


xt  =  xq  exp 


+ 


f  p(9s-)ds  +  [  cr(9s-)dBs+  f  f  G(z,9s-)ip(9s- ,dz,ds) 
JO  Jo  Jo  J\z\<l 

(6.1) 


■  f  f  H(z,  9s-),ip(9s- ,  dz,  ds) 

Jo  j\z\>l 

is  a  martingale  if  and  only  if 

M#t-)  +  L2(#t-)  +  /  [eG(z’et-)  -  G(z,9t-)  -  l}v(9t-,dz) 

+  [  -l}v(9t-,dz)=0,  M9t-  e  E. 

(2)  Furthermore  if  Let  in  (3.2)  is  replaced  by  M® 

dMQt  =  a(0t-  )dBt  +  [  G{z ,  0t-  )$(0t- ,  dz,  dt ) 

J\z\<  1 

-  [  H(z,0t-)^(Ot-,dz,dt), 

J\z  |>i 


(6.2) 


+ 


(6.3) 


then  the  solution  process  of  (3.1)  in  (3.18)  is  indeed  a  martingale  and  is  repre¬ 
sented  by 


xt  =  xq  exp 

pt 


\  [  cr2(9s-)ds+  [  cr(9s-)dBs 
z  Jo  Jo 

+  [  [  [ln(H(z,0s-)  +  1)  -  H(z,0s-)]v(6s-1dz)ds 

Jo  J\z\>l 
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+ 


+ 


+ 


/  /  ln(l  +  G(z,0s-))ip(8s-,dz,ds) 

Jo  J\z\<l 

[  [  [ln(l  +  G(z,6s-))  —  G(z,6s-)]v(6s- 

J 0  J\z\<l 

rr  ln(l  +  H(z,es-))$(Os-,dz,ds) 

Jo  J\z\>l 


,  dz)ds 


(6.4) 


(3)  Replacing  H(z ,  0S),  G(z,  Os)  and  L°s  in  (3.2)  by  eH^z,6s^  —  1,  eG^z,6s^  —  1  and  Me 
in  (6.3),  respectively,  the  solution  of  the  (IVP)  (3.1)  in  (3.20)  is  a  martingale  if 
and  only  if 


Xt  =  xq  exp 


+ 


[  a (0s-)dBa+  [  [  H(z,  6S- )v(6s- ,  dz)ds 
Jo  Jo  J\z\>l 

7;  [  cr2(0s-)ds+  f  f  [G(z,es~) +l]u(9s-,dz)ds 
^  J 0  Jo  J\z\<l 

[  [  G(z,0s-)ip(8s-,dz,ds)+  [  f  H(z,es-)tp(es-,dz,ds) 

Jo  J\z\<l  Jo  J\z\>l 


(6.5) 


(4)  V (t,  yt-  ,  6t- ,  xt- )  in  (4.13)  is  a  martingale  if  and  only  if  CV (t,  yt- ,  6t- ,  xt- )  is 
identically  equal  to  zero.  In  particular,  from  (4.34),  the  solution  process  (4.1) 
is  a  local  martingale  if  and  only  if  CV(t,  yt- ,  0t- ,  xt- )  in  (4.35)  is  identi¬ 
cally  zero  that  is  y(0t~)  +  H (z,  et- )v(9t- ,  dz)  +  EjeE\{et. }  KeZ  ~ 

l)\et_,j(yt-)]b{z\8t-,j)dz  =  0. 

We  introduce  and  recall  a  few  notations  necessary  for  presenting  the  next  lemma. 


Remark  6.1.  We  denote  $>t  a  positive  (P,  (Mt  VLt)te[0?T])-martingale  process  with 
initial  value  4>o  =  1.  In  fact  for  xo  =  1,  any  one  of  the  solution  processes  in 
Illustration  6.1  can  be  represented  by  4>t,  that  is  the  fundamental  solution  process 
of  linear  Levy- type  stochastic  differential  equations.  Moreover,  is  called  a  density 
process  of  a  probability  measure  P  with  respect  to  a  given  probability  measure  P. 


Based  on  a  Girsanov  theorem  for  Jacod  &  Shiryaev  (1987)  and  point  Bremaud 
(1981)  processes,  we  present  a  Girsanov- type  theorem  for  stochastic  hybrid  process 
described  by  (4.1).  We  highlight  the  effects  of  change  of  measures  on  both  time  and 
state  domains  of  decomposition  with  respect  to  (Lt)te[o,T]5  (Pn)n> o  and  (0t)te[o,Tp 
(Tn)  n>o  are  the  jump  times  in  Definition  2.1. 

Lemma  6.1  (Girsanov- type  Theorem).  Let  y  andY  be  piecewise  deterministic 
stochastic  processes  defined  on  [0 , T]  x  M  and  [0 ,T]  xRxl  into  M,  respectively. 
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£  =  z))mXm  is  a  M 171X171  -valued  and  Mt -predictable  process  defined  on  [0,T]  x 

M  into  M.  Let  us  consider  the  process  M[ 9  defined  by 


dM°=  f  Y  (e*-l  )N(dt,dz,{(Ot-,j)})+dL 

JzeR  Ar-  rr\  r n  \ 


(6.6) 


jeE\{et-} 


where  L\  is  defined  in  (3.2).  Furthermore ,  we  make  the  following  assumptions 

(HI)  £  =  (s,  ^))mxm:  {^i,j  (l/s))mxm  defined  in  (2.9),  tj  and  Y  satisfy  the  fol¬ 

lowing  conditions : 


<00  y  >  0 


(2,  s) Aij  (2/s-  )6(^  |  j)dzds 
eM 

/  rj(s,es-)/jL(03-)ds  <  00 
Jo 

f  [G(z,0s-) l\z\<!  Y  H(z,es-)l\z\>1][Y(Os-,z,s)  -  1  }v(Os-:dz)  <  00 
^^gm 

f  [Y(6s-,z,s)  -  l]2v{es-,dz)  <  00. 

<  Agm 


(6.7) 


(H2)  Ze£  6e  solution  process  of  the  following  linear  SDE , 


—  y+- 


V(t ,  0t)cr(6t)dBs  +  [  (Y (Ou  z ,  t)  -  l)^(0t,  dz,  dt) 

JzeR 


+  Y  f  -l}N{dt,dz,{{i,j)}) 

(i,j)eE2  JzeR 


Z0  =  1,  (6.8) 


where  N  =  N  —  7  defined  in  (4.14)  and  Zt  has  the  closed  form  representation 

Zt  =  exp  -[  Y(s,  9S- )2cr(9s-  )2ds  +  [  r](s,Os-)a(Os-)dBs 
.  Jo  z  Jo 

+  [  f  (Y(8s-,Z,S)  ~  1)^(0 s-,dz,ds) 

Jo  JzeR 

+  [  [  [\n(Y(6s-,z,s))  -  (Y(6s-,z,s) -l)]if(6s-,dz,ds) 

Jo  JzeR 

x  n  exp  /  0  - &j(s,z))b(z\i,j)\ij(ys-)dzds 
(i,j)eE2  '-Jo 

+  f  [  In (£itj(s,z))N(ds,dz,{(i,j)}) 

Jo  JzeR 


(6.9) 
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Therefore ,  from  Remark  6.1  and  under  a  local  equivalent  probability  measure  P  with 
density  process  Zt  with  respect  to  P,  the  following  hold : 

(1)  Bf  =  —  Jq  rj(s,  6s-  )cr(Os- )ds  +  Bt  is  a  Brownian  motion  for  each  6t-  G  E, 

(2)  up  (6t-  , )  =  Y (6t- ,  t,  )v(Qt- , )  P -almost  surely , 

(3)  p(dz,  {( i,j )})  =  £ij(t,z)b(z  \  i,j)\iij(yt)dz  P-almost  surely , 

(4)  Mf  defined  in  (6.6)  can  be  expressed  as  follows: 


dM l 


p(et-)  +  a2(et-)rj(t,et-)  +  [  G(z,et-)(Y(et~,z,t)-i)u(et-,dz) 


+  f  H (z,  9t-)Y (Ot- ,  z,  t)v(6t- ,  dz) 

+  E  [  [e'-lh  P(dz,{9t-,j}) 
jeE\{et-}JzeR 


dt  +  a(9t-)dBp 


+  [  G(z,0t-)jp(6t-,dz,dt)+  f  H(z,et-)^p(et-,dz,dt) 

+  T  [  [ez-l][N(dt,dz,  {Ot-,j})  ^p(dz,  {6t-,j})dt\. 
jeE\{et-}JzeR 


(6.10) 


Proof.  From  (6.8),  we  note  that  E[Zt  —  Zs  |  V  Ls]  =  0,  V  8,  t  G  [0,  T]  with  s  <  t. 
Hence,  Zt  is  a  local  martingale.  From  the  initial  condition  Zq  =  1  in  (6.8)  and 
Illustration  6.1,  Zt  is  a  density  process  of  P.  Moreover,  P(A)  =  fA  Zt(w)dP(w ),  for 
A  G  Ht  VLt.  Consequently,  P  is  a  local  equivalent  probability  measure  with  density 
Zt  relative  to  P.  From  the  definition  of  Bp  in  1,  it  is  obvious  that  it  is  a  Brownian 
motion  with  mean  rj(s ,  6S- )cr(0s- )ds  and  variance  t.  It  remains  to  show  that  Bp 
is  a  local  martingale  with  respect  to  P.  For  this  purpose,  we  use  (6.8)  and  apply 
Ito  formula  for  the  product  ZBP  Ladde  &  Ladde  (2013)  and  we  have, 

d(ZtBp)  =  ZtdBp  +  BpdZt  +  dZtdBp 
—  ZfdB  bJ^  dZt 
=  Zt[l  +  rj(t,et-)a(et-)Bp]dBt 


+  ZtBt 


I  ( 

J  z£R 


( Y(9t,z,t )  -  1  )ip(9t,dz,dt) 


+  Z  f  [£i,j(t,z)~l\N(dt,dz,{(i,j)}) 

U  A\^T?2  d Z£ R 


From  this,  we  conclude  that  Bp  is  a  P-continuous  local  martingale  with  quadr¬ 
atic  variation  t.  From  Levy  characterization  of  Brownian  motions,  Bp  is  a 
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P-standard  Brownian  motion.  This  establishes  1.  We  now  prove  that  vp(0t- ,  dz)  = 
Y (6t- ,  £,  z)v{0t-  ,  dz)  is  the  P-intensity  measure  of  ^(0t- ,  Knowing  that  P  and 
P  are  equivalent,  following  the  argument  Applebaum  (2009),  we  define  the  con¬ 
ditional  characteristic  function  for  the  Poisson  process  relative  to  the 

probability  measure  P  as  follows: 


Eb 


exp 


iu  /  ,dz,ds) 

Jo  JzeR 


l)]i/P(08,dz)  , 

(6.11) 


where  uP  is  an  intensity  measure  of  'ip  with  respect  to  P.  Using  the  closed  form 
expression  of  the  density  process  (6.9),  the  characteristic  function  in  (6.11)  is  also 
computed  as  follows: 


Eb 


exp 

=  E 

=  E 


iu 


nip(6s-,ds,dz) 
eM 

iu  /  ,dz,ds) 

Jo  J z CM 


exp 


Ht 


exp 

t 


[  [  (Y(03-,z,s)-l)i/>(03-,dz,ds) 

Jo  JzeR 

+  f  f  [ln(F(6>s- ,  z,  s))  —  Y(0S- ,  z,  s)  +  1  +  iu\ip(0s- ,  dz,  ds) 
Jo  JzeR 

[  [  (1  -Y(0s-,z,s))v(83-,dz)ds 

Jo  JzeR 

[  [  0'(Os-,z,s)~l)ip(Os-,dz,ds) 

Jo  J ze R 

+  f  pn(y  (Qs- ,  z,  s))  -  (Y (ds- ,  z,  s)  -  1)  +  iu]i/j(03,  dz,  ds) 
JzeR 

-  [  [  (Y {0S-  ,z,s)  -  1  )v(es- ,  dz)ds 
.  Jo  JzeR 


Mx 


=  exp 


x  E 


-  exp 


x  E 


exp 


n[ln(y  (0S- ,  z,  s))  +  iu\ip(0s- ,  dz ,  ds) 
gm 


Mx 


(6.12) 


We  note  that  fQ  J^eR[ln(F (0S-  ,  z,  s))  +  iu\ip(6s- ,  dz,  ds)  is  a  compound  Poisson  pro¬ 
cess.  From  Applebaum  (2009),  (6.12)  becomes 


Et 


exp 


=  exp 


iu  /  ^{0S-  ,ds,dz) 

Jo  Jz£R 

-[  [  (Y(0s-,z,s) -l)v(03-,dz)ds 

Jo  J  ze m 
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=  exp 


=  exp 


=  exp 


['  [  [eMY(0s,z,s))+iu  _  dz^ds 

0  J  z eM 

f  f  [Y(6S- ,  z,  s)eiu  -  Y(9S- ,  z,  s)\v(0s- ,  dz) 
Jo  J zeR 

[  [  Ke*“  -  l)\Y (' 9s - ,  z,  s)v{6a- ,  dz) 

Jo  JzeR 


(6.13) 


From  (6.11)  and  (6.13),  it  is  obvious  that  the  intensity  of  Levy  jump  poisson  measure 
with  respect  to  P,  is  v  =  Yu  P  —  almost  surely.  Based  on  the  proof  of  2,  the 
proof  of  3  can  be  reformulated,  analogously.  The  verification  of  (6.10)  follows  from 
algebraic  computations.  □ 


Remark  6.2.  We  recall  that  under  the  historical  probability  measure  P, 
( Pi,j)(i,j)eE 2  in  (2-5),  F(|i,  j)  in  Lemma  2.1  and  b(\i,j)  in  (4.3)  are  the  transition 
probability  matrix  of  the  embedded  Markov  chain,  the  sojourn  time  distribution  and 
the  log  jump  density,  respectively.  We  denote  (p^j)b,j)e^25  FP(| pj)  and  bp(\i,j ) 
the  transition  probability  matrix,  the  conditional  cumulative  distribution  of  sojourn 
times  and  the  density  of  the  log  of  jump  due  to  the  semi-Markov  process  from  state  i 
to  state  j  at  jump  time  Tn_i,  under  the  probability  measure  P.  Using  these  notions 
and  part  3  of  Lemma  6.1,  we  have 


bp(z  |  i,j)XPj(ys)  =  €ij(s,  z)b(z  \  with 

,p  (  x  p  fp(ys\i,j) 

\j(ys)  =  Pi,j- — ^ — PVT' 

1  -Ek^iPi,jFP(ys\hk) 


(6.14) 


We  further  remark  that  P  is  a  risk  neutral  measure,  if  the  process  Let  —  r(s)ds 
is  a  local  martingale  with  respect  to  P,  whenever  the  drift  coefficient  satisfies  the 
condition: 


-  r{t)  +  a2(dt-)ri{t,dt-)  +  [  G(z,dt-)(Y(dt-,z,t)  -  l)v(9t-  ,dz) 

J\z\<l 

+  f  H(z,9t-)Y(9t-,z,t)v{9t-,dz) 

+  [  ieZ -1hP(dz,{(°t-,j)})=  0.  (6.15) 

jeE\{et_} JzeR 


Given  the  2-variate  process  (rj(t,0t),Y(0t,  zj))  in  (6.8),  one  can  freely  choose  £. 
Hence,  for  each  choice  of  £,  one  gets  a  distinct  risk  neutral  measure.  Furthermore, 
by  the  application  of  the  first  and  the  second  fundamental  theorem  of  asset  pric¬ 
ing  Back  &  Pliska  (1991),  the  market  under  consideration  is  arbitrage  free  and 
incomplete. 
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Following  arguments  in  Elliott  et  al.  (2005),  Momeya  (2012)  and  Miyahara 
(1999),  we  define  two  particular  equivalent  martingale  measures,  namely  the  condi¬ 
tional  and  the  unconditional  minimum  entropy  martingale  measure,  respectively. 


6.1.  Conditional  minimum  entropy  martingale  measure 

We  define  the  conditional  minimum  entropy  martingale  measure  (CMEMM)  pricing 
Levy  and  semi-Markov  jump  risks.  In  the  absence  of  risk  associated  with  regime 
changes,  a  pricing  kernel  is  computed  through  a  random  Esscher  transform.  Without 
loss  in  generality,  we  assume  that  investors  always  know  past  and  future  market 
regimes.  Based  on  the  idea  in  Miyahara  et  al.  (2001),  we  define  the  process  Re  as 
follows: 

R?=  ft  E  [(ez-l  )N(ds,dz,{(0s-,j)})]+Ldt,  (6.16) 

Jo  Jz^jeE\{Os_} 

where  ,  n{t)  and  (3k  are  defined  in  (3.2),  (2.3)  and  (4.1),  respectively.  Picking  a 
locally  bounded  process  (&t)te[o,T],  we  modify  the  process  defined  in  (6.16)  as 

f  asdR°  =  f  [  E  as-[(ez -l)N(ds,dz,{(0s-,j)})] 

Jo  Jo  JzeR  jeE\{0s_} 


+  [  as-dLeq. 
Jo 


(6.17) 


In  the  following,  we  utilize  the  modified  process  (6.17)  to  formulate  a  dynamic 
process  for  the  asset  process. 


Definition  6.1.  Let  a  be  a  locally  bounded  process.  We  assume  that 
E[efo  as-dRds  |Ht]  <  oo,  V t  G  [0,T].  We  define  the  stochastic  processes  Za  and 
fc(s,  z,ds,dz)  as  follows: 


So<*s-dRes 


eJo 


E[efo  as-dRs  |  Mt\  ' 


Vt  >  0 


(6.18) 


and 


fc(s,  z,  ds,  dz) 


E  as-(ez  -  l)N(ds,dz,{(9s-,j)}), 

jeE\{es-} 


Vs  >  0,  z  G  M. 

(6.19) 


The  stochastic  process  defined  in  (6.18)  is  called  an  Esscher  transformation  with 
Esscher  parameters  (os)sG[0?t]- 


We  first  establish  preliminary  results  useful  for  finding  a  necessary  and  sufficient 
condition  under  which  the  probability  measure  Pa  with  density  relative  to  P  defined 
by  the  Esscher  transform  in  (6.18)  is  an  equivalent  martingale  measure  relative  to 
the  asset  price  process  xt  described  by  (4.1). 
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Lemma  6.2.  Under  Definition  4.1,  Remark  4.1  and  the  Esscher  parameter 
(as)sG[o5T]  in  Definition  6.1,  a  stochastic  process  exists  and  satisfies  the  following 
properties. 


(i) 


E[xf  |  Hr] 

n(t) 

=  n£ 

i= 0 


exp 


OLi 


rTi+ 1  r 

/  /  fc(s,  z,ds,dz) 

J  Pi  J  z  eM 


OiA 


i+ 1 


JT?+1  fi(s)ds 


(6.20) 


where 


fi(s)  =  a(i)n(di)  +  -a2(di)a(i) 


+ 


f  [e"WG(*A)  -  1  -  a(i)G(z,  9i)]v[9u  dz) 

J\z\<l 

f  i ]v(Qudz), 

J\Z\>1 


(6.21) 


for  i  e  1(0,  oo),  s  e  [0,  T}. 


(2) 


n 


T- 

exp  [a*  J^l+1  f  R  k(s ,  2:,  (is,  (P)]e^T* 


ao_  dMp 


(6.22) 


|  Ht]  ^  £?[exp[ai  J^+1  fzeR  k(s ,  z,  (is,  dz)]  |  6^  6i+1\ 
where  as-dM ^  =  as-dL ^  —  fi(s)ds. 

(3)  I  = 

(4)  Z?  =  zs  a  (P,  Ht  V  L)-ZocaZ  martingale. 

(5)  If  Pa  is  a  risk  neutral  measure  with  respect  to  Zff ,  then  under  Pa  we  have : 

(a)  Bpa  =  Bt  —  ft  as-a(ds-)ds,  is  a  Pa -standard  Brownian  motion  process. 

(b)  vpa(6s-,dz)  =  e^^,6,s-)1d^i>i)+G'(2:’0s-)1i^i<i]zy(6>s- ,  dz),  is  a  Pa -predict¬ 
able  compensator  of  the  Poisson  random  measure  ^(j, )  for  all  j  €  E. 

(c)  The  density  of  the  nth  jump  coefficient  f3n  is 

exp[P+1  fzeR  k(s’  z ’ ds’ 


E[exp[/r"+1  fzeR  k(s,  z,  ds,  dz)]  I  en,  en+i] 


Proof.  From  Definition  4.1,  0  =  To  <  T\  <  X2  <  *  *  *  <  Tn_i  are  the  regime  switch¬ 
ing  times  caused  by  the  semi-Markov  process  prior  to  t.  For  notational  convenience, 
we  denote  6L1  =  0q.  Under  the  assumption  of  the  lemma,  the  solution  process  of 
(4.1)  in  the  context  of  (4.2)  and  the  simple  return  process  (6.17)  exist  and  it  is 
represented  as 


x 


a 

t 


«(*)  r  ,T4+1  , 

T  exp  /  /  k(s,  z,  ds,  dz) 

,=f)  J  Ti  he  R 


f?+1 

eJTi 
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with  Po  =  xq  =  1.  For  t  G  [Tn,Tn+i],  from  the  independence  of  Levy  and  semi- 
Markov  processes,  we  have 


n—  1 

rTi+i  r 

E[x?  |Hr]  =  IJ 

E 

exp 

/  /  k(s,  z,ds,dz) 

2=0 

J  Ti  J  z  eiK. 

0i-i,0i 


xE[eJ 


s_dLiJ*as-dLes 


(6.23) 


This,  together  with  an  application  of  the  Levy  Kintchine  formula  0ksendal  &  Sulem 
(2005)  yields 


n{t) 

rE+ 1  r 

- 

1= 

II 

17 

exp 

/  /  k(s,  z,ds,dz) 

0i-iA 

2  =  0 

J  Ti  J  zeR 

x 


+ 


+ 


exp 


/: 


s-Kh)  +  2  a2(di)a2s- 


[  [e«s-G(^A)  _  1  -  as-G(z,9i)]v(6i,dz) 

J\z\<l 

L 


oas-  H(z,&i)  _ 


—  l]z s(6i,dz) 


This  completes  the  proof  of  (1).  For  the  proof  of  (2),  we  consider 


xt 

E[xa\nT]m 

From  (6.21)  and  (6.24),  we  obtain 


(6.24) 


E[xa  |Ht] 

_  nrlo  exP[/ff+1  LeR  z> dz)\  exP IJtT  ^s-dLT\ 

nrio  ^iexp  [Jt-+1  fzem  k(s’ z ’ ds’ dz )]  i  °j+ 1]  exp  [Jt-+1  fj(s)ds\ 

_  ttt  exp[/^i+1  fzeR  k(s,  z,  ds,  dz)}  exp[/^i+1  [as-dLesj  -  fi(s)ds]] 
i= 0  EiexP [JtT  fzeR  k(s’ Z ’  ds’  dzf  I  e3-l’0j\ 

(6.25) 

From  (6.3),  (6.17),  (6.21)  and  (6.25),  we  observe  that  as-dLesJ  —  fj(s)ds  has  a  form 
similar  to  (6.3),  that  is 

a3- dLesj  -  fj  ( s)ds  =  as-  dM°tj ,  (6.26) 
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with  coefficients  G  and  H  replaced  by  eG  —  1  and  eH  —  1,  respectively,  hence  estab¬ 
lishing  (2).  Using  (1),  (6.25)  and  (6.26),  we  further  remark  that 


E 


1 

^p 

E[x?\Mt] 

hAi'J' 

=  E 


pj  exP  [JtE  Lem  k(s’ z’ ds’ dz )]  exP  Ut-+1  as~dM P 


to  EiexP [fri+1  fze R  k(s’ z’  ds> dz )]  I  0i+i] 


yr  ^[exp [/^i+1  fz gR  k(s,  z,  ds,  dz)]  |  9j- 1, 9j] 

i=0  -E'[-^[exp[/y+1  f  R  k(s ,  z,  ds,  dz)]  |  0*,  0i+i]  |  0*_i,  0*] 


x  £ 


n(t) 


exp 


rT+i 

/,  [ 


'aa-dLej  -  fi(s)ds ] 


Mj1 


1 1  =  1,  for  *€[0,7], 


z=0 

which  establishes  (3).  For  the  proof  of  (4)  we  consider 

rTn 


E[xf  I IT] 


yr  exPK  ItE  fze K  k(S’  Z ’  ds’  dzf 

rg  OL  I  I  rj~i 

£[*?W]  i=„(s)+i  ^[exP[/Tii+1  Leu  k(s > ds> dz )]  I  °i  i> 


The  conditional  expectation  with  respect  to  Ht  V  Ls  yields 


E 


E[xf  |  H t] 


Mx1  V  Ls 


This  proves  (4).  Moreover,  from  (1),  (4)  and  (6.22),  Za  is  a  probability  density 
process  of  a  probability  measure  Pa  with  respect  to  P.  The  proof  of  statements  in 
(5(a))  and  (5(b))  of  (5)  follow  by  imitating  the  proofs  of  (1)  and  (2)  of  Lemma  6.1. 
We  only  establish  (5(c)).  For  B  C  and  t  G  [X/C,T/C+ 1).  In  fact 

Epa[lB]=E[lBzn 

=  E[E[lBZp  |  Hr]], 


=  E 


E 


exr 


VISE1  fzeRk(s,z,ds,dz)} 


n(t) 

n-  t 

to  ^[exp[/T;+1  fzeR  k(s,  z,  ds,  dz)]  \  0*_i,  0,] 
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=  E 


exP UtP  JzeR  fc(s>  dzj\ 

E[exp[J^+1  fzeR  k(s,  z,  ds,  dz )]  \9k-i,  Ok] 


n(t) 

x  n  e 

i=l,i^k 


exP Ut!+1  fzeR  fc(s’  ds ’  dzf 


E[exp[/ri+1  f  R  k(s,  z,  ds,  dz)} \9i-i,  9{ 


Oi-iA 


=  E 


1 B 


exP[irfcfc+1  Ler  fc(s> dzf\ 

E[exp[f^+1  fz£R  k(s,  z,  ds,  dz)}  |  9k- j,  6»fe]  _ 


exP[/T/+1 


1.  From  Radon 


Hence,  V5  6  B&,  T^fl#!  =  E\1b~  i.,1 

E[e*P[fT(+  /zGRfcO>Ms,cte)]  |  0fc-i,0fc] ' 

Nikodym  theorem  Jacod  &  Shiryaev  (1987),  the  density  of  (3k  under  P° 


exP [/t^  fzeR  h(s,z,ds,dz)] 

E[exp[ff*+1  fzeRk(s,z,ds,dz )]  |  0fc_i,0fc] 


.  This  completes  the  proof  of  the  lemma. 


is 


□ 


In  the  following  lemma,  we  provide  a  sufficient  condition  for  the  price  process  to 
be  a  (Pa,  (Ht  V  Lt)tG[0jT])-marf ingale.  The  result  obtained  will  be  used  to  derived 
the  martingale  condition  on  the  discounted  price  process. 


Lemma  6.3.  In  addition  to  assumptions  of  Lemma  6.2,  we  assume  that 
J^l >;L(i7(z,0s)  +  1  )ea^H(z'e^v(j,dZ)  <  oc,Vj  G  E.  Then  the  following  results 
hold : 


(1)  x  in  (4.1)  is  a  ( Pa ,  (H t  V  Lt)te[o ^-martingale  measure  provided  that 

n(9n)  +  ata* 1 2(9n)  +  [  G(z,  9n)[eatG<-z’en'1 * *  -  1  }v(0n,  dz) 

J\z\<l 

+  [  H{z,9n)eatH^e^p{9n,dz)  =  0,  (6'27) 

J\z\>\ 


Epa[(3n\en-i,en}  =  1,  Vi  e  (Tn,Tn+i),  in  e  1(0,  oo). 


(2)  The  discounted  price  process  Xt  =  e^o  Tsdsxt,  is  a  (Pa,(M.T  V  Lt)te[o,T])- 
martingale  if: 

p(Qn)  +  ata2(6n)  +  [  [G(z, 9n)eatG^z’6ni  -  G(z, 9n)}u(9n, dz) 

J\z\<l 

<  +f  [e^H{z,en)  _l]v{e^dz)  =  ru 

J\z\>l 

Epa  [[3n  |  9n-i,  9„\  =  1,  Vi  £  (Tn,  Tn+\),  VnS/(0,  oo). 

(3)  Let  a*  and  Pa  be  a  solution  process  of  equation  (6.27)  and  the  probability 

measure  associated  with  the  density  process  Za  ,  respectively.  Under  Pa  ,  the 
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process  Bfj.  in  (6.16)  could  be  expressed  as  follows: 

>e~[  “*-(*'-1  )N(dt,dz,{(6t-,j)})+rt-dt  +  a(6t-)dBp° 

JzeRjeE\{9t-} 


dR»  = 


+ 


+ 


f  G(z,dt-)[ip(9t-,dt,dz)  -  vpa  (i 9t-,dz)dt } 
J\z\<l 

f  h(z,  et- )  met- ,  dt,  dz)  -  vpa*  [et- ,  dz)dt] , 
J\z\>l 


with 


Epa*  [f3n  |  On- 1,  On\  —  1. 


Proof.  From  Radon  Nikodym  theorem,  xt  is  a  (Pa,  (H t  VLt)te[0,T])-marlingale  if 
and  only  if  xtZ^  is  a  (P,  (H t  V  Lt)tG[o,T])"marfingale-  From  (3.18)  and  (4.1) 


n{t) 

xtz ? = xsz«  n 

z=n(s)  +  l 


Aexp[/y+1  fzeRk(s,z,ds,dz)]eP  [“s  dMd+dLs^\ 


JTi  Jzem 

E[exp[J^+1  fzeWLk(s,z,ds,dz)}  \  6^-4 ,0*  \  0t- 1, 


with  Let  defined  as  follows: 


dLe„  = 


(6.28) 


ds 


p9s~)  -  ta2{0s-)  +  [  [ln(l  +  G(z,9s-)) -G(z,08-)]i/(08-,dz) 

Z  J \z\<l 

+  cr{Os-)Bs  +  f  ln(l  +  G(z,0a-))rp(63-,dz,ds) 

J\z\<l 

+  [  ln(l  +  H(z,  0s-))if(6s-  ,  dz,  ds).  (6.29) 

J \z  |>1 


'\z\>l 

From  (6.3),  (6.28)  and  (6.29),  we  have 
dA/fg*  H-  dPg 


P9S-) -^al-a2(9s-) -ta2(9s-)ds  -  f  [eas-H(z’es- )  -  l]v{6s- , dz) 
z  z  J\z\>i 

+  [  [as-G(z,9s-)-eas-G^ds-)  + 1 


+ in  (G(z,  9S- )  + 1)  -  G(z,  9S- )M  9S- ,  dz)ds 


+  {pLs-  +  \)<j(Os-)dBs 


+  [  [ois- G(z,  0S- )  +  ln(G(z,  6S- )  +  l)]ip(Os- ,  dz,  ds) 

J\z  l<i 

+  [  [as-H(z,0s-)  +ln(H(z,0s-)  +  1)]^0S~ ,  ch,  ds). 

Pd>i 


(6.30) 
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From  (3.18),  d[efo  d(as- Ms+Ll)j  =  e/o  with 


dL *  = 


+ 


M(^a-)  +  as-a2(ds~) 

[  [G(z,Os-)ea‘>-G<'z’e°-')  -  G(z,6s-)]v(0s-,dz)ds 

JUKI 


a.-H(z,et-) 


v(Qs-,dz) 


ds  -\-  <j{6s-  )(c^s-  ~\~  1  )dPs 


+ 


+ 


[  H(z,Os-)e 
J\z\>l 

f  l(G(z,Os-)  +  1  )ea.-G(*,e.-)  _  1  ]$(es_,dz,ds) 

J\z\<l 

f  [(H(z,es-)  +  I)e«s-H{z,es-)  _  1  }^es-,dz,ds). 
j\z  I>1 


(6.31) 


We  now  derive  conditions  under  which  xtZ^:  is  a  (P,  (H t  V  L£)te[0)T])-marfingale 
process.  XtZf  is  a  (P,  (Mt  V  L£)te[0;T])-marfhigale  process  if  and  only  if 

E[xtZt  |  Ht  V  Ls]  =xsZsa,  Vs, £  E  [0,T].  (6.32) 

Applying  Lemma  4.1  to  V(s,ys,6s,  Zsxs)  =  and  replacing  G,  id,  cr,  fi  and  fa 
by 

(G(z,9s)  +  l)eQsG7>0s)  —  1, 

(ff(zA)  +  l)e“*H(j!’®')  -  i, 

(as  +  1)<ts, 


M(^«)  +  oisa2{6s)  +  [  [G(z,  Os)easGG'e^  -  G(z,  6s)\v{6s,dz) 

J\z\<l 

-[  [easH{es) -l]v{6s,dz)  and 
Aexp[/J+1/2£R  fc(s,  2),  ds,  dz)] 


£[exp[/J+1  fzeR  k(s,  z,  ds,  dz )]  I  0i_i,  0»] 
respectively,  we  obtain 


xtZ?-xsZ“=  I  xu-Z« 


Js  Xu-K-  v^u-)  +Uu-a2{9u-) 

+  [  [G{z,  9u-)ea™- )  —  G(z,  6U- )]v{0u- ,  dz) 

JUKI 


7  i 

■'M>i 


|eau-H(*,flu_)  _  1]^(0„_  , 
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+  [  xu-ZZ-[(H(u,Ou-)  +  i)e«„- *(*.*„-)  _  i]v(pu_ 

J\z\>l 

+  [[  E  Xu~Zu-(eZ  ~l)N{duidzdeu-,j}) 

Js  JzER  jeE\{6u_y 


dz) 


du 


+  sum  of  martingale  terms  . 


Taking  the  conditional  expectation,  we  obtain 


E[xtZ, “  -  |  HT  V  La] 

=  f  E[xu-Z%-  |HrVLs]  +  a„-cr2(6»„-) 

J  S 

+  [  [G(z,6u-)eav-G(-z’ev-'>  -  G(z,9u-)]v(9u-,dz) 

-  f  H(z,0u-)ea«-H^-K(0u-,dz) 

ft  V  au-xu-Z^(ez -l)N(du,dz,{(6u-,j)}) 

J  S  J  ZGM  rr,\  rn  ~\ 


du 


+  E 


jeE\{eu-} 


Mx1  V  Ls 


=  0,  Vs, £  G  [0, T], 


(6.33) 


for  any  s,  £  and  for  small  As  s,  £  =  s  +  As  G  (Tn,  Tn+i)  for  some  n  G  7(1,  oo).  This 
together  with  (6.33)  yields 

H(9S)  +  ascr 2(6s)  +  [  [G(z,  0s)easG('z’es'>  -  G(z,  9s)\v(0s,dz) 

JUKI 


+  [  H(s,Os)ea°H(z’e^v(0s,dz)Aa  =  O,  Vs  e  (Tn,  Tn+1),  n  e  7(1,  oo). 
J|Z|>1 

(6.34) 


Lastly,  we  assume  [s,  £]  =  [Tn,  Tn+i].  When  As  is  small.  There  is  one  regime  change 
[s,£]  at  £  =  Tn.  Using  (6.34)  and  (6.33)  becomes 


E 


rTn+i  n 

/  /  E  «„-[xu-Z“_]|HTVLr„ 

JTn  j  zER  Ae-  rp_\  rn  1 


ZTnxnE 


jeE\{Ou _} 
^ean(/3n-i) 


,L[ead/^n  i)  |6>n_!,^ 


'  |^n— 1 5  |I®T  V  LTri 


=  0,  Vn  G  7(0,  oo) 

—  1  =  0,  Vn  G  7(0,  oo) 
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E 


E 


faeOnVn-V 


E[ean(Pn  1)  I  en„i,6n 


9n  — 1 >  9n 


—  1  =  0,  VnG  /( 0,  oo), 


Pn  exp[/T”+1  fzeR  k(s,  z,  ds,  dz )] 


.  E[exP[fTf+1  fzeR  fc(s’ ds’ dz )  I  ^n]] 


0n-i,0n  -1  =  0,  Vn  G  7(0, oo), 


Ep<x  [pn  |  $«_i ,  $n]  1  =  0,  VnG/( 0,  oo). 

(6.35) 

This  completes  the  proof  of  (1).  (2)  is  a  direct  consequence  of  (1)  whenever  fJt{0s-) 
is  replaced  by  n(6s-)  —  rs.  For  the  proof  of  (3),  we  use  (6.27)  to  derive  the  risk 
neutral  dynamic  of  the  process  R°  defined  in  (6.16).  We  denote  Bp  and  vp  the 
standard  Brownian  motion  and  the  intensity  process  of  the  Poisson  process  ijj  under 
the  probability  measure  Pa  ,  respectively.  From  Lemma  6.2  (5(a)),  solving  for  B 
in  Bp  =  Bt  —  f0  a*_cr(6s-)ds ,  adding  and  subtracting  vp  inside  the  Poisson 
integrals,  we  obtain: 

dRt  =  [  E  (e*  -  l)N(ds,dz,  {(0t-,j)})  +  n(0t-)dt  +  a(6t-)dBt 

J ze E  m  r/>  1 


jeE\{et-} 


+  [  G(z,9t-)'tp(0t- ,dz,dt)  +  f  H(z,et-)i)(6t-,dz,dt) 

J\z\<\  J\z\>\ 

=  [  E  (ez-l  )N(dt,dz,m-,j)}) 

J ZeR  •  _  m  !o  X 


jeE\{et-} 


n{9t- )  +  v2  iPt-  )at- 


+ 


+ 


f  G(z,et-)[upo:  (9t- ,  dz)  —  v(9t- ,  dz)] 
J\z\<l 

[  H(z,9t-)ppa  ( 9t-,dz ) 

J\z\>l 


dt  +  a  [9t-  )dBp 

+  f  G{z,9t-)[ip(9t-,dt,dz)  -  ppa  ( 9t-,dz)dt] 

J\z\<l 

+  [  H(z,9t)[ip(9t,dt,dz)  -  ppa  (9t,dz)dt],  Vi  G  [Tn,Tn+1].  (6.36) 

J\Z\>1 

From  Lemma  6.2  (5(a)),  one  gets  vp“  (j,dz )  =  e^H^z’^1(\:‘i>1')+GGd)1\z\ <dv{^j^dz). 
Hence, 

,d  ~  [  E  (e*-l )N(dt,dz,{(6t-,j)})  + 

JzeRjeE\{et- } 


dRt  = 


/j,(9t-)dt  +  a2(9t-)a*- 


+ 


[  G(z,9t-)[ea°t-G(z’et-]  -  1  }v(9t-,dz)dt  +  a(9t-)dBI 

J\Z\<1 
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+  /  H{z16t-)eXet-H^z'dt  \(6t-  ,dz)dt 

+  [  G(z ,  et- )  ,  dt ,  dz)  -  z/**  (0t- ,  cfe)dt] 

+  [  H(z,Ot)[ip(Ot,dt,  dz)  -  vpa  (9t,dz)dt],  (6.37) 

J\z\>l 

where  a*  satisfies  the  condition  (1).  Therefore,  (6.37)  becomes 
dRt  =  [  E  (^-1  )N(ds,dz,{(0t-,j)})  +  a(0t-)dBpa* 

JzeRjeE\{et _} 

+  f  G(z16t-)[ip(0t-  ,dt,dz)  —  vp  (0t~,dz)dt\ 

J\z\<l 

+  [  H(z,6t-)[fi(Ot-,dfidz)-upa\ot-,dz)dt}1  VtG  [Tn,Tn+1], 

(6.38) 

with 

[/?n  |  0n-lj  0n]  =  1, 

which  proves  (3).  This  establishes  the  lemma.  □ 


In  the  next  remark,  we  introduce  a  particular  case  of  Rf  corresponding  to  the 
simple  return  process  Miyahara  et  al.  (2001)  and  we  present  a  few  properties  of 
conditional  entropies  Fujiwara  (2003). 


Remark  6.3.  Let  Pi  and  P2  be  two  absolutely  continuous  probability  measures 
relative  to  P.  We  recall  three  important  properties  of  conditional  entropies  Miyahara 
(1999) 

(1)  74-VLt(Pi|P)>0. 

(2)  w^(Pi  I  P)  <  I  P),  if  K  C  eT  V  Lt. 

(3)  If  Pi  is  a  (P,  (H t  V  Lt)te[o,T])"absolutely  continuous  martingale  measure,  and 

P2  is  a  probability  measure  equivalent  to  P  such  that  ln(^-)  is  integrable  with 
respect  to  Pi,  then  I  P)  >  EPl  [In (dP2  |  dP)  |  Hr]. 

We  now  state  and  prove  the  conditional  minimum  entropy  property  of  the  mar¬ 
tingale  measure  Pa  when  Rf  is  the  simple  return  process  of  xt  in  Remark  6.3. 


Lemma  6.4.  Let  Q  G  M(x,P)  =  {Q  <C  P  :  x  is  a  (Q,  (Ht  V 
Local  martingale} .  Let  Pa  be  defined  as  in  Definition  6.1  with  a*  solution  pro¬ 
cess  of  (6.28).  Then  the  following  inequality  holds : 


0/H71 

'PffiyVLy 


(Q\P)>n 


My1 

My  VlLy 


(P“*|P), 


VQe  M(x,P), 


where  Wgv Ht 


zs  £/ie  conditional  entropy  defined  in  Elliott  et  al.  (2005). 
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Proof.  We  prove  the  lemma  in  two  steps.  The  first  step  consists  in  minimizing  the 
conditional  relative  entropy  of  any  probability  measure  Q  in  the  set  M(x,P).  From 
(6.18)  and  Remark  6.3,  one  notes  that  Z can  also  be  expressed  as  follows: 


Z?  = 


E[e 


Jo  *s-dR0s 


|Ht]‘ 


(6.39) 


By  definition  of  an  absolutely  continuous  local  martingale  measure,  the  discounted 
stock  price  x  =  e~ fo  rs~dsxt  is  a  (Q,Wt  V  £t)te[o5T] -local  martingale  process.  The 
simple  return  processes  Rt  and  Rt  are  associated  with  the  price  process  xt  and 
the  discounted  price  process  xt  with  respect  to  (4.1),  defined  by  dRd  =  and 

dR?  =  4?^-,  respectively.  Q  £  J\A(x,P)  implies  that  x  is  a  (P,  Mp  V  LA-  mart  ingale. 
Hence,  R%  is  a  Q-local  martingale  process.  Furthermore,  fQ  as-dR°s  is  a  local  Q- 
martingale.  As  a  Q-local  martingale,  a*_  dR°  is  therefore  integrable  with  respect 
to  Q.  From  (6.39)  and  (6.21),  we  have 


where 


In (Zf)=  f*  a*-  dRg  —  f  g(s)ds, 
Jo  Jo 


g(t)  =  In (E[eti<-dk°  \  Hr]). 


(6.40) 


(6.41) 


From  (6.40),  ln(Z"  )  is  integrable  with  respect  to  Q  as  a  sum  of  two  integrable  terms. 
Let  (tn)ne/( o,oo)  be  a  local  sequence  of  increasing  stopping  times  with  limn^00  tn  = 
T,  associated  with  the  local  martingale  as-dRes.  By  definition  of  local  sequences, 
the  process  J^nM  a s=dRd  is  a  Q-martingale.  Hence,  for  any  t  £  [0,T],  we  have 


KU  JQ  I P)  >  KWz,,,„  (Q I  n  (Remark  6-3) 


(  dP^ 

\ 

>  Eq 

!  "e 

ElrVLtAtn  / 

Mx1 

,  (Remark  6.3).  (6.42) 


From  (6.40),  we  have 


f  dPa 

\ 

Eq 

In 

y  dP 

Ht  VLtAtn  / 

Mt 

=  EQ[ln(Z?;tJ\UT} 


=  EQ 
=  Eq 


ots-  dRs  |  Mj1 


+  pQ  [d(t  A  tn)  |  Ut] 


asdRs  |  Ht  V  Lq 


+  Pq  [git  A  tn)  |  H t] 


=  Eq  [g(t  A  tn)  |Mt], 
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since  J^Mn  a s-dR(s)  is  a  (H t  V  Lt)  t e [o , r] -mart ingale .  We  note  that,  Eq  [g ( t  A 
in)|HT]|  <  Eq[\g(T)\  |  HT]  =  g(T)yt  e  [0,T],  since  from  (6.41),  g  is  HT- 
measurable.  Hence,  by  the  Dominated  Convergence  theorem,  we  have 


lim  Eq  [g(t  A  tn)  |  HT]  =  Eq  [g(T)  |  HT]  =  \g(T) | .  (6.43) 

n— >  oo 

Taking  the  limit  in  (6.42),  we  obtain 

KlvlrW  I  P)  >  EQ[g{T)  |  Hr]  =  g(T).  (6.44) 

The  second  step  of  the  proof  consists  in  showing  that  the  conditional  relative  entropy 
of  the  random  Esscher  transform  achieves  the  minimum  value  in  (6.44).  Using  (6.39) 
and  the  Pa  -martingale  property  of  Rt  the  relative  entropy  of  Pa  is  computed  as 
follows: 

Ht 


<UT(pa*\  p)  =  Epa 


In 


dPa 

dP 


=  Epa * 


as-dRs  |  Hr  V  Lq 


+  Epa*  [g(T)  |  Hr],  (from  (6.39)) 

=  Epa*  [g(T)  |  WT]  =  g(T),  (. Ret  is  a  Pa  martingale). 


From  (6.44),  the  lemma  follows. 


□ 


6.2.  Unconditional  minimum  entropy  martingale  measure 

We  will  define  an  equivalent  martingale  probability  measure  and  we  will  estab¬ 
lish  that  it  has  the  unconditional  minimum  entropy  martingale  measure  property. 
(A ij(t))mxm  is  the  intensity  matrix  of  the  semi-Markov  process  0  from  (2.9)  and 
N  is  the  point  process  defined  in  Definition  4.4. 

Definition  6.2.  Let  Q  be  a  local  absolutely  continuous  probability  measure  with 
respect  to  the  historical  probability  measure  P  on  the  filtered  measurable  space 
(U,  Ht  VLt,  (Ht  VLt)tG[o,T]);  Q  and  P,  denote  the  regular  versions  of  the  conditional 
probabilites  P(|Ht)  and  Q(|Ht)  over  H t  V  Lt-  Zf  denotes  a  local  martingale 
process  with  initial  value  1,  representing  the  density  process  of  Q  with  respect  to 
P.  {£i,j{t))mxm  denotes  a  matrix  entries  with  predictable  processes.  Moreover,  the 
rows  add  up  to  0  and  satisfy  Xq*  j)eE2  Jo  (s)^,j(s)l^s  <  oo.  N  and  7  are  the 
processes  defined  in  Definition  4.4. 

We  first  recall  a  decomposition  theorem  Momeya  (2012)  and  we  establish  a 
Girsanov-type  lemma  necessary  in  the  proof  of  the  UMEMM  property. 

Lemma  6.5.  Let  Q ,  P,  Z®,  Q(w, .),  (£i,j)mxmj  AT,  N  and  P(w , .)  be  proce¬ 
sses  and  probability  measures  defined  in  Definitions  6.2  and  4.4.  The  following 
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claims  hold : 

(1)  There  exist  two  density  processes  ZtL  and  Z fjd  such  that 


z ?  =  Zt  x  ZT 


H 


with 


dQ 

dP 


=  ZtL  and 

t  dP 


(6.45) 


=  Zr 


H 
T  * 


Mx 


Lt  VHt 

(2)  If  (A ij(t))mXm  is  the  matrix  with  conditional  intensity  of  the  semi-Markov 
process  0  in  (2.9)  and  A ij(t)  ^  0 fit  G  [0 ,  T],  then  the  following  claims  are 
equivalent : 

(a) 

dQ 


dP 


=  ,  where  ZM  solves  the  SDE: 


H 


Ht 

H 


dZ?  =  Zt 


H  E 

(■ i,j)eE 2 


-1 


A  i,j{t  ) 


N(dt,R,{(iJ)}),  Zq  =  1. 


(b)  Under  probability  measure  Q,  with  density  process  Zp ,  the  point  process  M 
has  conditional  intensities  matrix  (€ij(t))mxm- 


Proof.  The  proof  of  (1)  follows  closely  Momeya  (2012).  As  for  (2),  we  note  that 
(2(a))  =>  (2(b))  follows  from  the  proof  of  Lemma  6.1.  We  now  aim  at  proving 
that  (2(b))  =>  (2(a)).  From  Definition  6.2  and  (2(b)),  Zp  and  N(,  R,  {(i,  j)})  — 
y(R,  {(i,  j)})  are  (P,  (H£)te[0?T])-marfingale  processes.  From  the  martingale  rep¬ 
resentation  property  of  N(t,R,{(i,j)})  =  N(t,R,{(i,j)})  -  j(R,{(iJ)})  = 
N(t,R,  {(i,j)})  —  A there  exists  an  m  x  m  matrix  of  Ht-predictable  processes 
(slJ)mxm  such  that 

dZtH= 

(ij)eE* 

As  Zp  >  0  P-almost  surely,  there  exists  an  m  x  m  matrix  of  predictable  processes 
sf3  satisfying  s\,J  =  Zpsf3.  Hence, 

dZ?  =  Z»  Y,  %iN(dt,R,{{i,j)}). 

( hj)eE ^ 

From  Lemma  6.1,  the  matrix  of  conditional  Q-intensities  of  JV(,  R,  {(i,  j)})  is 
\ij(t)(l  Psf3).  One  needs  to  prove  that  the  conditional  intensity  of  7V(,R,  {(i,j)}) 
with  respect  to  Q  is  £ij(t),\/i,j  G  1(1,  m).  Hence,  equating  both  matrices  and 
solving  for  s M  yields 

+  Slp)  =  €i,j(t),  \/t  G  [0, T\,  ( i,j )  G  E2,  and  hence, 


s+  = 


-1 


^ i,j  (t)  _ 


Vt  G  [0, T],  (iJ)eE2 
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Therefore,  ZH  is  solution  process  of  the  SDE 


dZtH 


=ztH  y 

(ij)eE*  L 


-1  + 


KA1  ). 


[N(dt,R,{(hj)})  ~  tijil  )dt\  ^o  =  l- 


From  Lemma  6.1,  the  intensity  matrix  of  the  semi-Markov  process  0  under  the 
probability  measure  Q  is  This  completes  the  proof  of  (2)  and  hence 

the  lemma.  □ 


We  define  a  density  process  which  we  prove  is  the  unconditional  minimum 
entropy  martingale  measure. 


Definition  6.3.  Let  Pa  ^  be  a  risk  neutral  measure  with  density  Zt  =  Z™  x 


Ti 


where  is  introduced  in  Definition  6.1,  with  Re  the  simple  return  process  of 
defined  in  (6.16).  a*  is  the  solution  process  of  (6.28)  and  is  solution  of  the  SDE 


dzf  =  z\_  £ 


t-  Z.MeEP&dM'p. 


dP 


So  a*  dRs 


eJ  o 


ItVLt  E(edo  "s-dfls|ElT) 


exp 


t)  ( i,j)eE 2 


[  ))Vi(s  )ds 

Jo 


■  Tin  (^(s-))N(ds,R,{(i,j)}) 
Jo 


We  also  define  a  functional  F  as  follows: 


mij))  =  e 


9(E)  +  E 

(■ iJ)eE 2  L 


[  (1  £i,j  (s  )ds 

JO 


[T\n  (6,i(5-))JV(d5,R,{(i,j)}) 

Jo 

[  (1  —  £i,j(s  (s  )ds 

Jo 

I1  ln(ei,i(5-))JV(d5,R,{(i,j)}) 

Jo 


II  exp 


where  g  is  defined  in  (6.41). 


We  will  next  show  that  under  a  particular  choice  of  £,  Pa  ^  has  the  uncondi¬ 
tional  minimum  entropy  martingale  measure  property. 

Lemma  6.6.  We  denote  Pa  ^  and  F  the  risk  neutral  measure  and  the  functional 
from  Definition  6.3,  respectively.  If  Q  is  a  (Ht  VLt)te[o,T]  risk  neutral  measure  and 
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( CtJ)mxm  minimizes  the  functional  F,  then  the  following  holds: 


mxm 


I P)  >  nMTWLT(P^^  I P). 


Proof.  Let  Q  be  a  risk  neutral  measure.  By  definition  of  risk  neutral  measures,  Q 
is  locally  absolutely  continuous  with  respect  to  P.  From  Lemma  6.5,  there  exists  a 
process  Zff  and  a  process  Z^1  such  that  ^p|HTVLt  =  % t  x  -  From  Lemma  6.5, 
we  have 


Zf=  P  exp 

rT 

/  (1  —  Ci,j(s  {Vs-  )ds  E 

/  In(e(a-))JV(ds,R,{(t,j)}) 

(*,i)es2 

Jo 

Jo 

for  some  mxm  matrix  of  H^-predictable  processes  as  in  Definition  6.2. 


=  E[Z*fZ* (In (Z£))  +  ZlfZ*  In (Zf )] 


=  E\E\ZJfZ^  In (Z£)  |  Ht]  +  E[Z^Zfi  In (Zf )  |  HT]] 
=  E[Z* E[ZJf  ln(Z^)  |  Ht]  +  Zf  In (Zf )] 


=  E[Z»H™IvLt(Q  |  p)  +  z»  ln(Z")] 
>  P[Z#^v£t(P“  |  P)  +  Z*  In (Zf )] 


—  E  p(p)  +  ^  '  /  (1  £i,j(s  ))^i,j(ys~)ds 

U  Pc  7^2  ^0 


x 


=  n&j) 

>  (definition  of  £) 


(p“*-«ip), 


which  proves  the  result. 


□ 
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6.3.  Siu  and  Yang  Kernel  pricing  all  risks 

Let  a  be  a  piecewise  constant  stochastic  process.  We  define  the  density  process  Z f 
of  a  probability  measure,  Pa,  on  the  filtration  (Mt  V  L t)te[o,T p  with  an  Esscher 
transform  with  parameter  a.  The  pricing  kernel  discussed  here  is  based  on  the  work 
in  Siu  &  Yang  (2009),  in  the  context  of  a  Markov  switching  asset  price  process. 


Definition  6.4.  Let  Za  be  the  following  stochastic  process: 


dPa 

~dP 


Hit  VLt 


=  Z? 


E 

1 


e~fo  a(s)dRes 

E[e-f?a--dRi  \90,yo\ 


Ht  V  L t 


if  Vt  G  (0,  T] 
if  t  =  0, 


where  Ret  is  the  log  price  process  in  (6.16),  induced  by  x  from  (4.1)  in  the  context 
of  solution  (4.2). 


We  note  that  from  Lemma  5.1,  one  can  retrieve  any  particular  scalar  con¬ 
ditional  characteristic  function  from  the  vector  characteristic  function  as  fol¬ 
lows:  T(iq  £,  y,  j,  x)  =  exp[irdn(x)](exp(/J+?/  M(u,  s)ds)  •  ej:  1),  where  e#t  = 
(l(9t=i,  l0t=2, . . . ,  l0t=m)T.  at  is  the  Esscher  parameter  process  associated  with  the 
probability  measure  Pa. 


Lemma  6.7.  Let  Zj*  be  the  process  in  Definition  6.4.  is  an  almost  surely  positive 
martingale  with  unitary  expectation. 


Proof.  We  first  prove  that  Z*f  is  a  martingale.  Let  0  <  s  <  t 


E[z ; 


t  I  ““S 


VLJ  =  E 


=  E 


E 


e-tf«s-dR°s 


E[e~fo  as-dK\y0je0,L0} 

So  as-dRs 


It  V  L t 


VL 


E[e~fo  as-dRes\y0,d0,L0\ 


V  L, 


(Hs  V  Ls  c  Ht  V  Lt) 


=  Zf. 


Therefore,  Zf  is  a  martingale.  It  follows  that  Zf  has  unitary  expectation, 
E(Zf)  =  E[E(Zf)  |  H0  V  L0]  =  E(Zf)  =  1. 


Noting  that  Zf  is  an  almost  surely  positive  process  by  construction,  the  lemma 
follows.  □ 


From  the  preceding  lemma,  Zf  is  a  density  process.  Hence,  the  Esscher  transform 
in  (6.4)  defines  a  probability  measure  Pa  equivalent  to  P.  It  remains  to  show  that 
Pa  is  a  martingale  measure  under  a  certain  condition  specified  in  the  next  lemma. 

Lemma  6.8.  Let  Zf  be  from  Definition  6.4  and  x  as  defined  in  (4.1).  M(u,y) 
and  M(u,y)  are  defined  in  (5.4)  with  modified  log  price  process  defined  by 


1550052-52 


Option  Pricing  with  a  Levy-Type  Stochastic  Dynamic  Model 


dRf  =  OLt-  ln(/?n(t))  +  at-dLet  and  dRet  =  at-  ln(/3n(t))  -  rt-dt  +  (at-  +  l)dL?» 
respectively ,  with  Le  defined  in  (3.2). 

Xt  =  e_^o  rs-dsxt  is  a  (Pa,  (HI*  V  ^t)te[o,T]) -martingale  if  and  only  if 

r  p+y 


exp 


•>  1 


■  rt-ey 

/  M(—i,s)d, 

J  v 

'  rt+y 

/  M(—i,s)ds 
Jv 


■e*0,l  )  =0  Vt  G  [0, T],  G  E.  (6.46) 


Proof.  Let  0  <  u  <  t.  From  Momeya  (2012),  Siu  &  Yang  (2009)  and  by  the  abstract 
Bayes  rule  (Jacod  &  Shiryaev  1987),  we  have 

E(Z?e~t »‘r-ds^|HuVL„) 


=  e~ttr'-daxt 


(6.47) 

(6.48) 


E[e~rrs-ds+K(a(s  )+! )dRds  |H„VL„] 

1  £[e/»  |H„VL„]  ' 

Hence,  e_^o  r»~dsXt  is  a  (_P“,  (EH/  V  f-n)te[o,T})- martingale  if  and  only  if 

^-r^-+1)-?l^V^]=l  Vu,  f  G  [0,T]. 

From  Lemma  5.1  applied  to  =  at-  ln(/3n(t))+at-dLj  and  =  at-  ln(/3„(t))  — 
rt-dt  +  (at-  +  l)dLf  .  respectively  and  on  account  of  the  Markov  property  and  the 
homogeneity  of  the  process  (0,y),  the  numerator  and  the  denominator  of  (6.49) 
becomes 


(6.49) 


E[e 


-f*r,-d,+f*[a(,~)+l]dLt 


VL„ 


v 


=  (  exp 


and 


/  ry+t-u  _  \ 

\i  M(-i,s)dsj  ■  ee0,'. 

E[et o*"(s  )dL?  |H„  VLJ  =  ^exp^y  M(-i,  s)ds'j  ■  eg0, 1^, 


respectively,  where, 

-^dp,q{  L  £/) 


~r(q)  +  (oiq  +  1  )fi(q)  +  \{otq  +  l)2cr2(<?)  +  f  [e(a<!+|)G(z’9)  -  l]i/(g,  dz) 

+  [  [e<'olq+1')G<'z’q) -l-(aq  +  l)G(z,q)]v(q,dz)  +  \q,q(y)  if p  =  q 

J\z  I<1 

\q,p(y)  /  e(aQ+1)^5(2;  |  q^p^dz  otherwise 
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and 


Mp,q(  'f'll/) 

aqy{q)  +  \°?(j2(q)  +  f  [ e aiG(z^  -  1  -  aqG(z,  q)]v(q ,  dz) 

1  J\z\<  1 


+  f  [ea«G^  -  1  ]i/(g,  dz)  +  X q,q(y)  if  p  =  g 
•'N>i 

A q,p(y)  /  eaqZb(z\  q,p)dz  otherwise. 


Hence,  (6.49)  becomes 

ry+t-u 


exp 


&0o  i  1 


ny-\-t-u 

/  M(—i,s)ds 

Jy 

ry+t-U 

/  M(—i,s)ds 
Jy 


e90,l)—Q,  Vrt,tG[0,T],  \/ 6U  e  E , 

(6.50) 


which  completes  the  proof  of  the  lemma. 


□ 


7.  Option  Pricing  Formulas 

In  this  section,  we  price  a  European  style  call  option  within  the  risk  neutral  pricing 
theory  Schachermayer  (2010).  We  denote  Q  an  equivalent  martingale  measure  of 
the  historical  probability  measure  P,  relative  to  the  price  process  x  in  (4.1).  We 
derive  a  PIDE  extending  the  PDE  in  Black  &  Scholes  (1973)  satisfied  by  European 
call  prices.  We  also  describe  how  two  existing  pricing  methods  blend  seamlessly  in 
the  context  of  this  paper. 

Definition  7.1.  Let  S  be  a  function  in  P2(f7,Q)  defined  on  M+  x  M+  into 
M  representing  the  payoff  of  a  contingent  claim;  Q  is  a  risk  neutral  probabil¬ 
ity  measure  of  the  price  process  x  defined  by  (4.1)  with  respect  to  the  histori¬ 
cal  probability  measure  P;  if  is  a  non- negative  real  number  denoting  the  strike 
price  of  a  European  type  option  contract  with  maturity  T;  xt  denotes  the  asset 
price  value  at  maturity;  C  is  the  Q-risk  neutral  option  price  function  defined  on 
[0,T]  x  M+  x  M+  x  [0,T]  x  E  x  M+  into  M+  and  V  denotes  the  discounted  option 
price  process  defined  by  V  ( t ,  T,  K,  yt,  6t,  xt)  =  e~fo  rs-dsC(t,  T,  K,  yt,  Qt,  xt). 

Lemma  7.1.  Let  S  be  a  random  variable  representing  the  payoff  of  a  general  Euro¬ 
pean  style  contingent  claim  with  maturity  T  and  strike  price  K  in  Definition  7.1; 
let  Q  be  the  risk  neutral  measure  defined  in  Definition  7.1  and  C  is  the  Q-risk  neu¬ 
tral  option  price  of  a  contingent  claim.  Then ,  the  Q-risk  neutral  option  price  C  of 
a  European  contingent  claim  with  maturity  T,  strike  price  K  and  payoff  S  can  be 
expressed 

C(t,T,K,yt,0t,Xt)=EQ(e-Pr‘dsS(xT,K)\yt,et,xt).  (7.1) 
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Proof.  From  Schachermayer  (2010),  the  Q-risk  neutral  option  price  C  at  time  t  is 
given  by 

C(t,  T,  K,yt,8t,xt)  =  Eq  [e~T  r‘dsS(xT,  K)  |  Ht  V  Lt], 

We  note  from  Lemma  2.3  that  the  triplet  (y,0,x)  is  Markovian,  hence 
C(t,  T,  K,  yt,  9t,  Xt)  =  EQ[e~P r°dsS(xT,  K)  \  yu  9t,  xt], 
which  proves  the  result.  □ 


A  partial  integro  differential  equation  (PIDE)  satisfied  by  a  European  style 
contingent  claim  with  maturity  T  and  payoff  H  is  presented  in  the  next  lemma. 


Lemma  7.2.  Let  Q ,  C  and  V  be  the  risk  neutral  measure ,  the  Q-risk  neutral  option 
price  function  and  the  discounted  option  price  process  defined  in  7.1,  respectively. 
Then  V  satisfies  the  following  system  of  PIDE 


dV  dV 

~d7  +  ~d^  +  Kj)Xs~ 


dV  1  2...  2 

+  r {,)x- 


d^V 

dx2 


L 


V(s,T,K,ys,6s,xs-  +xs-G(z,j )) 


N<i  L 

dV 

E •,  K ,  ys ,  j .  xs— )  G(z,  j)xs— 


v{j,  dz ) 


+  f  [V(s,T,K,ys,9s,xs-  +  xs-H(z,0s)) 

J\z\>l 

-  V (s,  T,  K,  ys,  0s,xs- )}v(6s,dz) 

+  /  E,hj(.ys)V(s,T,K,ys,j,xs-ez)b{z\i,j)dz 
JzeR 


-t~  T ,  K ,  ys—  i  i,  xs— ) A (y)  0, 

with  terminal  condition 

V(T,  T,  K,  yT,  0T  =  j ,  xT)  =  e~£  r^dsS(xT,  K),  for  j  e  E. 


Proof.  From  (7.1),  the  discounted  price  process  could  be  expressed  as  follows: 

V(t,  T,  K,  y,  j,  x)  =  e~fo  r^dsC(t,  T,  K,  y,  j,  x) 

=  E?{e-h*.n  r^dsS(xT,  K)  \  yu  9t,  xt)  (7.2) 

V  is  a  (Q,Lt  V  Hb))-Martingale  since  it  is  a  Q-conditional  expectation.  We  use  the 
law  of  iterated  expectation  and  u  <  t  to  prove  it  as  follows: 

E{V (t,  T,  K,  y,  j,  x)  |  L„  V  H„) 

=  E(e~fo  r^dsC(t,  T,  K,  y,  j,  x)  \  Lu  V 
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=  E[e~K  r^^dsE(e~  P  r^dsS(x,  K)  \  Lt  V  Ht)  |  Lu  V  H„] 

=  E[E(e~  P  r^dsS(x,  K)\Lt\J  Ht)  |  L„  V  H„] 

=  E[e-P  r^dsS(x,  K)\LUW  H»] 

=  eSo  r{es)dsE^e-  fj  r(e.)dsg^  ^)  |  £„  V  Hu] 

=  V(t,T,  K,y,j,x). 

From  Ito  differential  formula  in  Lemma  4.1,  we  have 

dV (£,  T,  K ,  yt,Ot,  x(t))  =  AV (£,  T,  if,  yt- ,  0t- ,  x(t~))dt  +  Martingale  Terms  . 


As  V  is  a  martingale,  the  first  term  vanishes  and  (4.14),  the  following  PIDE  is 
obtained: 

8V  OV  ,  8V  1  2/n  N  2d2V 

aT  +  %  +  )x-~  a£  +  2"  )a!-  a? 


+ 


s,T,  K,ys,6s,xs- 


+  xs- G(z ,  0S))  -  V" (s,  T,  if,  2/s,  a;s- ) 


-  G(z,  9s)xs- 


dV' 

dx 


v(Qs,dz) 


+  [  [V(s,T,  K,ys,es,xs-  +xs-H(z,9s))  -  V(s,T,K,ys,6s,xs-)]v(6s,dz) 

J\z\>l 

+  /  y2^ij(ys)V(s,T,K,ys,j,xs-ez)b(z\i,j)dz 
Jzem  jzjLi 

~V(s,T,K,ys-,i,xs-)X$s_j(ys)  =  0,  Vi  e  [0,T]. 

Hence,  the  proof  is  complete.  □ 


Definition  7.2.  Let  C  be  a  continuous  function  on  [0,T]  x  M+  x  M  x  [0,T]  x  E  x 
M+  into  M+  representing  the  modified  European  call  option  price;  let  T  be  the 
characteristic  function  of  C  with  respect  to  its  third  variable,  and  k  denotes  the 
logarithm  of  the  positive  real  number  if  in  Definition  7.1. 

Remark  7.1.  Assuming  a  deterministic  interest  rate  r,  a  closed  form  formula  for 
the  Fourier  transform  of  a  modified  vanilla  European  call  option  price  is  known 
Carr  &  Madan  (1999).  Let  us  denote  C,  C  and  y  the  European  call  price,  the 
modified  European  call  price  of  Carr  and  Madan  type  and  a  positive  real  number, 
respectively,  for  the  payoff  function  of  a  European  call  option  S(xt,K)  =  {xt  — 
if)+  =  (eln(XT)  —  efc)+,  with  k  =  ln(if).  Further  assume  that 

oo 

\C(t,T,  k,  j,y)\dk  <  oo,  Vj  e 
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Then  the  modified  European  call  price  defined  by  Carr  and  Madan  is  expressed  as 
follows: 

C(t,  T,  k,  yu  eu  x(t))  =  evEQ(e~  S?  r(t)d.^(xT)  _  ek)+  |  ^  ^  (7.3) 


In  the  following  lemma,  we  recall  the  characteristic  function  of  the  modified 
European  Carr  and  Madan  type  call  option  (Carr  &  Madan  1999). 


Lemma  7.3.  If  T  and  ^(u,t,yo,6o,x(0))  are  the  Fourier  transform  of  C  defined 
in  Definition  7.2  and  the  conditional  characteristic  function  of  the  log  price  process 
defined  in  Lemma  5.1,  respectively ,  then  we  have 


r(T,w,yo,0o  =  j,x0) 

e~fo  rs-ds 


(rj  +  iw)(  1  +  77  +  iw) 


wk  -  ix0(  1  +  rj),  T,  7/0, 6>0  =  j,  £0).  (7-4) 


Proof.  The  proof  can  be  found  in  Carr  &  Madan  (1999). 


□ 


Remark  7.2.  We  note  that  in  the  case  of  a  regime  switching  interest  rate  Momeya 
(2012)  uses  the  Carr  and  Madan  type  transformation  to  obtain  the  characteristic 
function  of  European  call  option  prices.  The  formula  in  Carr  &  Madan  (1999)  is 
based  on  the  characteristic  function  of  occupation  times  which  is  known  in  closed 
form  when  market  states  are  described  by  a  Markov  Chain.  We  have  derived  the 
characteristic  function  of  the  occupation  times  in  Corollary  5.1  which  allows  us  to 
extend  the  results  in  Elliott  &  Osakwe  (2006)  and  Momeya  (2012)  when  market 
states  are  described  by  a  semi-Markov  process. 


In  the  context  of  a  price  process  driven  by  the  Brownian  motion  (Ghosh  & 
Goswami  2009),  an  integral  option  price  formula  is  obtained.  In  the  following  result, 
we  present  a  similar  pricing  formula  (Ghosh  &  Goswami  2009)  in  the  context  of 
(4.1),  where  we  assume  that  /|  is  the  density  of  the  increment  of  the  log  price 
process  in  an  interval  of  length  s,  whenever  the  semi-Markov  process  is  in  state  j 
for  any  j  E  7(1,  m)  =  E. 


Lemma  7.4.  An  integral  option  pricing  formula  in  the  context  of  model  (4.1)  is 
represented  by  the  following  formula : 


C(t,  yu  9t,xt)  =  P(t,  Vu  0t)C9t  (t,  T,  K,  xt) 


■  Q(t,  yt,  St)  [  < 

Jo 


er^up(t,yt,9t) 


pOO 

/  C(t  +  u,0,j,xt)ck 
Jo 


dx , 
(7.5) 

with  P(t,yt,6t )  =  ’  1  ~  PdiVt,0t)  =  Q(t,yt,0t),  p(t,yt,6t )  = 

/i-Ffa7|ef)t)  and  +  u,yt,j,xt)  =  T,et+u=j,j^0t-  cd  +  u,yt,j,xt)P{\n(x/St)), 
where  F(-  \  0t~)  and  /(•  |  Qt~)  are  defined  in  Remark  2.1.  CBt  is  the  Black-Scholes 
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option  price  when  the  market  is  in  state  Ot  and  C(t ,  ytlOt,xt)  is  short  hand  notation 

forC(t,T,K,yt,Ot,xt). 

Proof.  The  lemma  follows  by  imitating  the  proof  of  Theorem  3.1  of  Ghosh  & 
Goswami  (2009).  Let  V(t)  =  v(t^ytl0t,xt)  defined  as  in  Lemma  4.1,  using  the  risk 
neutral  pricing  formula,  the  tower  Law  of  expectations,  the  identity  1  =  lTn(t)+1<i  + 
lTn(t)+i>i  and  the  notations 

E[V(t)\s,ya,0a,x3]  =Es[V(t)], 

E[V(t)  |  s,ys,9s,xs,Tn{t)+1  <T}=  EfT[V{t)\,t 
E[V(t)\s,ys,es,xs,Tn{t)+1  >T]=EfT[V(t)] 

E\V  {t)  I  u ,  y$,  Uw)  9  s,  9u,  xs, •i'u\  ES,U\V (t)] 

E[V(t)  |  s,u,ys,yu,0s,6u,xs,xu,Tn(t)+1  >  T]  =  E>pV(t)} 

E[V(t)\s,s  +  u,ys,ys+u,9s,9s+u,xs,xs+u,Tn(s)  =  yt  +  u]  =  ETs=T[V{t)\, 

V  s,  t,  u  £  M+  and  s  <  t,  we  obtain 
C (L  Uti  @t->  Xt) 

=  Et-  [eP  r{-e°)ds(xT  -  K)+] 

=  Et-  [E[eP  r^ds(xT  -  K)+  |  yu  9t,  xt ,  Tn(t)+1]] 

=  Et-[l{TnW+l)>TE[eSP^-)ds{xT  _  K)+  \yt,et,xt,TnW+1]] 

+  E[lTn(tH1<TE[ePrVs-)ds(xT  -K)+\yt,et,xt,Tn(t)+1}] 

=  P(Tn(t)+ i  >  T  |  yt- ,  9t- ,  xt-  )Epr[E[C(t ,  yt- ,  9t- ,  xt- )  |  Pn(t)+ 1]] 

+  P((Tn^t)+i  <  T)  |  yt- ,  9t- ,  )£'t-T[£'(C(t,  yt- ,  9t- ,  xt- )  |  Tn(t)+i]] 

=  P(Tn(t)+ 1  >  P I  J/t-  >  1  )-®[C(t,  t/t- ,  0(- ,  £t- )  |  Tn(t)+i  >  T] 

+  P{Pn(t)  + 1  <  ^  |  Vt-  >  ®t~  !  )-®[C(t,  2/t~  )  @t~  >  xt~  )  I  ^n(i)  +  l  <  r] 

=  P(Tn(t)+ 1  >  P I  J/t- )  ®t~ !  2-'*-  )CSt  (t,  Xt)  +  P(Tn^t)+i  >  T  |  yt- ,  ,  a:t- ) 

x  £'t-_r[£'(C(t, yt- ,  9t-  ,xt-)\  rn(t)  =  yt  +  u}] 

=  P(t,  yt,  9t)C6t  (t,  xt)  +  Q(t,  yt,9  +  t) 

[T~t 

x  /  =  yt  +  u\du 

Jo 

=  P(t,  Vt,  9t)Ce*  (t,  xt)  +  Q(t,  yt,  9t)  [  p(t,  yt,  9t)Et+u,t 

Jo 

x  =yt+u]\ 
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=  P(t,yudt)Ce*(t,xt)  +  Q(t,yt,9t)  [  *  Et[er^up(T  -  u,yt,  9t) 

Jo 

x  El=yt+u[ef^r{^-)ds(xT  -  K)+]]du 

=  P(t,yt,Ot)Cet(t,xt)  +  Q(t,yt,9t)  f  er(0t)>(T  -  u,yt,Ot)du 

Jo 

poo 

x  T]  E[C(t  +  u,yt+u  =  0,Ot+u,xt-)\Tn(t)=yt  +  u] 

JO  n 


Ot  +  u=j 

j^ot 


x  fl(\n(x/St))dudx 

P(t,yt,Ot)Cet(t,xt)  +Q(t,yt,6t)  [  er(-et)up{T  -  u,  yt,  0t) 

Jo 

pOO 

x  /  C(t  +  7i,  0,  j,  xt)dudx. 

Jo 


□ 


In  the  following  section,  we  provide  a  numerical  application  to  exhibit  the  role, 
scope  and  significance  of  the  developed  results.  Furthermore,  based  on  the  options 
on  DJIA  quoted  on  March  3,  2008,  we  also  present  calibration  results.  The  data  set 
was  retrieved  from  Deville  (2007). 


8.  Numerical  Applications 

We  present  simulated  European  style  call  option  prices  and  implied  volatilities 
derived  from  option  prices  generated  by  a  semi-Markov  switching  exponential  Levy 
pocess. 


8.1.  Simulations 

Option  prices  derived  in  our  framework  are  more  flexible  than  option  prices  derived 
from  the  more  commonly  used  Markov  regime  switching  exponential  Levy  pro¬ 
cesses.  We  also  demonstrate  that  prices  simulated  in  our  framework  produce  implied 
volatilities  consistent  with  stylized  facts  of  the  option  market,  namely,  the  smile  and 
the  smirk.  The  risk  neutral  measure  considered  in  this  section  is  the  conditional 
minimum  entropy  martingale  measure  developed  in  Sec.  6.  Simulations  of  option 
prices  and  implied  volatilities  will  be  performed  from  semi-Markov  regime  switch¬ 
ing  Black-Scholes  (SMBS),  semi-Markov  regime  switching  Merton  jump  diffusion 
(SMMJD)  and  semi-Markov  regime  switching  normal  inverse  Gaussian  (SMNIG) 
processes.  The  choice  of  these  three  Levy  processes  is  rooted  in  the  fact  that  they 
highlight  three  important  classes  of  Levy  processes:  Diffusion  processes,  jump  diffu¬ 
sion  processes  and  processes  with  infinite  activity.  At  each  market  state  we  assume 
j,  H(j ,  z)  =  G(j,  z)  =  z  in  (4.2)  and  the  SMBS  has  no  Levy  jump  component,  hence 
its  Levy  measure  is  z/(j,  •)  =  0,  Vj  G  E.  The  diffusion  coefficient  at  each  state  j  is 
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a(j),  V  j  G  E.  At  each  state  j,  the  SMMJD  has  Levy  jumps  arriving  at  a  Poisson 
rate  with  mean  r](j);  its  jumps  size  is  normally  distributed  with  mean  and  standard 
deviation  p(j)  and  S(j),  respectively  for  any  state  j  G  E.  From  Tankov  (2003),  it’s 
Levy  measure  is 


v{j,  dz) 


Ail 

5(j)V2n 


exp 


1  (z-p(j))2 

2  8*(j) 


dz. 


At  each  state  j  G  E ,  the  SMNIG  is  a  time  changed  Browian  motion  with  drift  rate 
/x(j).  The  new  time  scale  is  defined  by  an  inverse  Gaussian  process  with  parameters 
S(j)  and  p2(j)  —  The  SMNIG  is  a  pure  jump  process  and  consequently  does 

not  have  diffusion  component  ie  a(j)  =  0,  V  j  G  E.  When  the  market  is  in  state 
j  G  E,  the  characteristic  function  of  the  SMNIG,  is  expressed  as  follows: 


0(j,  u)  =  exp  [Sy/p2(j)  -  (pj)  +  m)2  -  y/p2(j)  -  rf{j)],  with  i  =  y/^1. 


Under  the  risk  neutral  measure,  the  price  jumps  at  times  of  regime  changes  Kij  are 
assumed  to  follow  the  distributions: 


1  - 
1 


^  1  +  £ij 


with  probability  pij 

with  probability  1  —  pij  —  qij 

with  probability  q^j. 


(8.1) 


The  risk  neutral  condition  (6.28)  requires  that  Epa*  (/3n  \  6n-i,  0n)  =  1  ,Vn  G 
/( l,oo).  This  implies  that 


Pi,j=Qi,j ,  V(i,j)e£'2. 


(8.2) 


From  (8.1)  and  (8.2),  the  characteristic  function  of  In (ftij)  reads 

E[eiu ’“(**•*>]  =  bij(e_i“ln(1+e<’j)  +  eiu(1-ln£i’i}  -  2)  +  1], 


£i,j  e  [0)1])  PiJ  € 


°G 


,  V(i,j)£E2 


(8.3) 


where  the  coefficient  of  u  is  i  =  and  the  subscript  “i”  of  ft,  e  and  p  stand 

for  the  ith  state  of  the  semi-Markov  process.  The  sojourn  time  distribution  of 
the  semi-Markov  process  is  assumed  to  be  a  piecewise  exponential  distribution 
approximating  a  two-parameter  Weibull  distribution.  It  is  worth  noting  that  piece- 
wise  exponential  distributions  are  known  to  have  piecewise  constant  intensities. 
Knowing  that  Riemann  integrable  functions  can  be  approximated  almost  every¬ 
where  by  step  functions,  we  can  simulate  prices  in  our  framework  with  a  wide 
range  of  sojourn  times.  Let  T*  be  the  time  horizon  of  the  market,  with  T  <  T*, 
where  T  is  any  option  contract  maturity  time.  We  consider  the  following  partition 
0  =  ao  <  a\  <  •  •  •  <  clm- i  =  T*  of  [0,  T*]  and  we  denote  clm  =  oo.  From  the  def¬ 
inition  of  the  Weibull  distribution  and  (2.9),  the  Weibull  intensity  function  can  be 
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approximated  by  the  following  piecewise  constant  semi-Markov  intensity  function 
associated  with  the  partition  (a)^Q1 


\j(Vs)  —  < 


PiJ  53  71  (  -T  )  1  [ak,ak+1)(Vs)  if  *  ±  j 

k=0  ^  ^  ^  ' 


^i,j{ys) 

j=l,j^i 


otherwise 


M—l 


=  < 


k= 0 


“i,i  53  WO*  ll\au,ak+1){ys) 

otherwise, 


(8.4) 


^i,j(ys) 

j=l,j^i 


Vs  E  [0,T],  with  =Pi,j^  and  a£  =  °fe+°fc+1 ,  Vfc  E  1(0,  M  —  2).  Figure  2  shows 
that  for  the  particular  parameters  chosen,  the  piecewise  constant  approximation  of 
the  Weibull  intensity  is  quite  close  to  the  actual  Weibull  intensity.  In  a  market  with 
Markov  regimes,  intensities  are  assumed  to  be  constant.  This  is  quite  a  restrictive 
assumption,  since  it  implies  that  as  time  goes  by,  the  market  propensity  for  chang¬ 
ing  regimes  stays  constant  even  when  market  conditions  change.  It  is  indeed  more 
realistic  to  allow  the  intensity  matrix  to  be  time  dependent  and  vary  accross  the 
lifespan  of  the  derivative  contract.  This  is  precisely  the  advantage  of  semi-Markov 
over  Markov  market  regimes,  as  semi-Markov  intensities  depend  on  the  time  spent 
on  the  current  state.  It  is  for  instance  well  known  that  markets  pre  and  post  2007 
financial  crisis  are  quite  different  as  the  volatility  has  gone  up  a  notch  and  prices  have 
become  more  erratic  Manda  (2010).  We  therefore  expect  switching  rates  between 
market  regimes  in  the  post  crisis  market  to  be  different.  Option  price  simulations  are 
performed  through  inversion  of  the  Fourier  transform  of  the  modified  option  price  in 
Carr  &  Madan  (1999).  Inversion  is  performed  via  the  fractional  Fourier  transform 
described  in  Chourdakis  (2004),  as  it  allows  a  greater  flexibility  in  the  choice  of 
the  log  strike  grid  parameter  and  the  frequency  grid  parameter.  The  size  of  the  log 
strike  grid  and  the  frequency  grid  is  set  at  3V  =  4096.  The  infinite  integral  in  the 
frequency  domain  is  truncated  at  512  and  the  domain  of  ln(^)  is  set  to  be  [—1,  +1]. 
Under  this  setting,  the  fractional  parameter  is  e  =  ^  x  x  «  9.72  x  10-6. 
Throughout  this  section,  we  assume  that  the  market  has  m  =  4  distinct  regimes, 
xt  =  $100,  the  interest  rate  is  constant  across  all  states  r  =  (0.05,0.05,0.05,0.05), 
the  initial  backward  recurrence  time  y  —  0.1  year,  the  time  horizon  of  the  mar¬ 
ket  is  T*  =  5  years,  the  time  to  maturity  is  set  to  T  =  0.5  year  for  option  price 
simulations  and  T  G  [0, 1.5]  for  “Implied  volatility”  surface  simulations.  The  param¬ 
eters  specific  to  our  framework,  not  accounted  for  in  markets  with  Markov  regimes 
are  the  backward  recurrence  time  y  at  time  of  pricing  and  the  shape  parameter 
vector  ($i)i£E  of  the  sojourn  time  distributions.  It  is  therefore  of  interest  to  see 
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Fig.  1.  Weibull  intensity  versus  sojourn  time  distribution  intensity. 


Fig.  2.  Semi-Markov  regime  switching  NIG  option  prices  at  all  market  regimes  for  different  shape 
parameters  values. 


how  much  difference  those  parameters  induce  when  option  prices  from  a  market  with 
semi-Markov  regimes  are  compared  with  option  prices  from  a  market  with  Markov 
regimes  ceteris  paribus.  The  partition  associated  with  the  piecewise  constant  semi- 
Markov  conditional  intensity  considered  (8.4)  is  a  regular  grid  of  [0,  T*]  where  M  = 
61  with  a/e  =  k  ]^_1 ,  k  E  1(0,  M  —  1).  Figure  1  shows  how  well  A  approximates 
weibull  intensities.  This  choice  of  M  is  rooted  in  the  fact  that  financial  performances 
of  firms  are  usually  made  available  on  a  quarterly  basis.  We  therefore  consider  three 
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Semi  Markov  switching  M.TD  with  <  1. 


Markov  switching  M.TD  with  i),  =  1. 


Strikes 


Semi  Markov  switching  MTD  with  0t  >  1. 


Fig.  3.  Semi-Markov  MJD  option  prices  at  all  market  regimes  for  different  shape  parameters  values. 


0  100  200  300  400  500  600  0  100  200  300  400  500  600 


Semi  Maikov  vs  Maikov  NIG:  Regime  3  Semi  Markov  vs  Markov  NIG:  Regime  4 


Fig.  4.  Semi-Markov  NIG  option  prices  at  all  market  regimes  for  different  shape  parameters  values. 


periods  within  each  quarter:  The  earning  month  where  the  market  likely  (over) reacts 
to  reports,  the  post  earning  month  with  a  possible  market  correction  and  the  pre 
earning  month  where  the  market  is  preparing  for  the  next  earning  report.  The 
intensity  of  the  semi-Markov  process  6t  of  Definition  2.1  is  assumed  constant  within 
each  month.  Parameter  values  chosen  for  simulations  of  option  prices  and  volatility 
surfaces  are  summarized  in  Tables  1-4,  respectively.  Option  price  simulations  from 
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Fig.  5.  Comparison  of  simulated  European  call  option  prices  derived  from  Markov  and  semi-Markov 
switching  MJD  models. 


Markov  regime  switching  models  and  semi-Markov  switching  model  with  increasing 
intensities  are  performed  using  ==»  (1, 1, 1, 1)  and  ^  =  (2,  3,  8,  5),  respectively.  In 
Figs.  2  and  3,  one  notices  that  deep  in-the-money  option  prices  for  ^  <  1  ,  Vi  G 
1(1,  m)  and  y  =  0.1  year  flirt  with  the  $120  mark,  while  prices  induced  by  $*  >  1 
and  =  1,  Vi  G  1(1,  m)  lie  below  $100.  In  the  presented  simulated  prices,  it  appears 
that  the  initial  backward  recurrence  time  and  the  shape  parameter  of  sojourn  time 
distribution  do  induce  noticeable  price  differences  between  options  in  Markov  and 
semi-Markov  market  states.  We  also  note  price  differences  between  market  states 
as  shown  in  Figs.  2  and  3.  On  the  other  hand,  we  note  that  by  changing  the  shape 
parameter  of  the  sojourn  time  distribution,  semi-Markov  prices  prove  to  be  more 
exible  than  Markov  prices  as  shown  in  Figs.  4  and  5.  One  recalls  that  the  higher 
the  shape  parameter  values,  the  smaller  the  mean  sojourn  time  of  the  market  at 
each  regime  and  the  smaller  the  regime  change  risk.  This  partially  explains  why 
deep  in-the-money  options  in  Figs.  4  and  5  are  cheaper  as  $  gets  larger.  As  far  as 
implied  volatilities  are  concerned,  smiles,  smirks  and  other  shapes  are  seamlessly 
reproduced  as  evidenced  in  Fig.  6. 

8.2.  Calibration 

The  Dow  Jones  Industrial  Average  Index  (DJIA)  is  a  performance  indicator  of 
30  large  public  US  companies  that  accounts  for  approximately  20%  of  the  U.S. 
market.  We  consider  options  quotes  on  the  Dow  Jones  Industrial  Average  Index 
(DJX)  on  March  3,  2008,  as  the  global  financial  crisis  was  in  full  effect.  We  calibrate 
the  parameters  of  our  models  (SMBS  and  Markov  BS)  to  market  option  prices  in 
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Semi-Markov  BS  Implied  Volatility  Surface  Semi- Maikov  MID  Implied  Volatility  Surface  Semi- Maikov  NIG  Implied  Volatility  Surface 


Fig.  6.  Implied  volatility  surfaces  from  option  prices  generated  by  semi-Markov  switching  BS,  MJD 
and  NIG  respectively  at  the  first  market  state. 


Table  5.  We  use  residual  sum  of  squares  (RSS)  to  compare  the  fit  of  both  models. 
The  model  parameters  sought  are  chosen  so  as  to  minimize  the  sum  of  squared 
differences  between  market  options  and  model  option  prices.  We  assume  that  m  —  2; 
the  option  contract  has  maturity  time  T  =  47  days;  the  underlying  stock  does  not 
pay  dividends  and  is  currently  valued  at  S  =  $122;  the  sojourn  time  intensities  are 
assumed  to  be  a  piecewise  constant  approximation  of  Weibull  intensities  as  defined 
in  (8.4);  the  fractional  parameter  is  given  in  Sec.  8.1  and  the  current  backward 
recurrence  time  is  assumed  to  be  y  =  252  days.  Calibration  results  obtained  are 


Table  5.  DJX  option  quotes  of  March  3,  2008. 


Strike  price 

Strike  S/K 

Prices  (T  =  47) 

98 

1.25 

24.43 

99 

1.23 

23.40 

100 

1.22 

22.50 

101 

1.21 

21.55 

102 

1.20 

20.63 

103 

1.19 

19.68 

104 

1.18 

18.75 

105 

1.16 

17.83 

106 

1.15 

16.90 

107 

1.14 

15.98 

108 

1.13 

15.10 

109 

1.12 

14.23 

110 

1.11 

13.33 

111 

1.10 

12.45 

112 

1.09 

11.63 

113 

1.08 

10.78 
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Table  5.  ( Continued ) 


Strike  price 

Strike  S/K 

Prices  (T  =  47) 

114 

1.07 

9.95 

115 

1.06 

9.18 

116 

1.05 

8.40 

117 

1.04 

7.68 

118 

1.04 

6.93 

119 

1.03 

6.23 

120 

1.02 

5.58 

121 

1.01 

4.95 

122 

1.00 

4.35 

123 

0.99 

3.80 

124 

0.99 

3.25 

125 

0.98 

2.74 

126 

0.97 

2.28 

127 

0.96 

1.90 

128 

0.95 

1.52 

Table  6.  Markov  BS  and  SMBS  parameter  estimates. 


Model 

(do,  di) 

(#0,#l) 

(eo,ei) 

(po,Pi) 

(do,  di) 

RSS 

SMBS 

(70.85,  1.523) 

(26,  0.202) 

(0.008,  0.55) 

(0.001,  0.02) 

(0.212,  0.18) 

2.07 

Markov  BS 

(79.03,  9.363) 

(1.  1) 

(0.176,  0.175) 

3.73 

Regime  1:  Markov  BS  prices  vs  Observed  prices  Regime  2:  Markov  BS  prices  vs  Observed  prices 


Strikes  Strikes 

Fig.  7.  Markov  BS:  Calibrated  versus  observed  prices. 


summarized  in  Table  6  and  plotted  in  Figs.  7  and  8.  One  note  that  the  RSS  is 
considerably  lower  for  semi-Markov  regime  switching  models.  In  fact,  the  RSS  is 
reduced  by  44%.  The  fit  of  SMBS  is  visibly  better  as  shown  in  Figs.  7  and  8. 
Moreover,  the  scale  parameters  of  the  Weibull  shows  that  market  option  prices  are 
best  explained  by  semi-Markov  market  regimes  with  =  26,  $2  =  0.202.  These 
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Regime  1:  SMBS  calibrated  vs  Observed  prices  Regime  2:  SMBS  calibrated  vs  Observed  prices 


Strikes  Strikes 

Fig.  8.  SMBS:  Calibrated  versus  observed  prices. 


shape  parameter  values  show  that  state  1  and  2  have  increasing  and  decreasing 
intensities.  More  precisely,  the  longer  the  market  stays  in  state  1,  the  more  likely  it 
is  to  switch  to  state  2  and  the  longer  the  market  stays  in  state  2  the  more  likely  it 
is  to  stay  in  state  2,  respectively.  A  look  at  the  discrete  jump  size  also  shows  quite 
a  significant  risk  neutral  probability  (2%)  of  drop  or  jump  in  price  of  more  than 
50%  when  the  market  is  in  state  2.  This  is  in  line  with  the  widespread  fear  and  the 
observed  erratic  market  prices  during  the  financial  crisis.  Indeed,  on  September  29, 
2008,  the  DJIA  suffered  its  biggest  one-day  price  drop  (7%)  and  a  few  month  later, 
November  13,  2008,  the  market  registered  its  biggest  one-day  price  jump  (11%). 
More  importantly,  the  Dow  Jones  suffered  a  whopping  33%  drop  in  2008. 


9.  Conclusion 

In  this  work,  we  proposed  a  theoretical  framework  for  asset  price  models  by  devel¬ 
oping  a  closed  form  an  exponential  Levy  under  a  semi-Markov  switching  process.  In 
addition,  we  established  a  few  distributional  properties  of  the  price  process  in  (4.1). 
We  derived  an  infinitesimal  generator  of  the  triplet  (yt,0t,xt)  and  the  character¬ 
istic  function  of  the  log  price  process  Let  which  allowed  us  to  get  a  closed  form 
expression  for  the  characteristic  function  of  plain  vanilla  European  style  call  option 
prices.  Such  a  characteristic  function  could  prove  to  be  useful  in  the  calibration  of 
model  parameters  to  market  options  prices  through  inverse  Fourier  transform,  as 
demonstrated  in  Sec.  8.  Two  problems  of  interest  are  currently  being  investigated: 
(1)  an  evaluation  of  the  effects  of  the  backward  recurrence  time,  the  sojourn  time 
distribution  and  the  jumps  distribution  on  option  prices  induced  by  semi-Markov 
switching  asset  price  models;  and  (2)  the  pricing  of  exotic  options  in  semi-Markov 
market  regimes. 
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